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A. INTRODUCTION 


A rocket in flight has three linear and three angular degrees of freedom. A 
mathematical model m be devised which describes the motions within the six 
degrees of freedom ana this '6-D* model may be programmed for solution by a 
digital computer. Given a set of initial conditions and tables of rocket characteris- 
tics, the equations which make up the model may be evaluated repeatedly with 
respect to time to give an accurate simulation of &e flight of a rocket. 

The ’6-D' model, however, has two deficiencies: 

1. Significant effort is required in preparing the input data to be used by 
the computer program. A specimen '6-D' simulation shows that over 
700 ii^ut data items are required. When trajectories must be predicted 
for a new type of rocket, only part of this data set may be available 
from the manufacturer. The remaining items must be calculated fiv,m 
aerodynamic theory. All the data must then be transcribed in the exact 
form reqiiired. As a result, preparation of the ii^)Ut data usually takes 
at least two man-weeks. 

2. Each '6-D' computer run may require as much as four hours of expensive 
computer time. 

For these reasons, models utilizing various sub-sets of the six degrees of freedom 
are used in trajectory simulation. A '3-D' model with only linear degrees of 
freedom is especially attractive, since the coefficients for the angular degrees 
of freedom are the most difficult to determine and the angular equations are the 
most time consuming for the computer to evaluate. 

Of course, the '3-D' model is less accurate than the '6-D' model. This is because 
the model lacks angular motions and the thrust vector orientation is assiimed to 
be aligned with the velocity vector. Unless the angle of attack is zero, this is not 
true. 


Figure 1 shows a typical angle of attack versus time history of an unguided rocket. 
The time scale is divided into tlireu periods. In Period I the rocket has a finite 
an^e of attack and is said to be imtrimmed. In Period n the angle of attack has 
been trimmed out to zero and in Period in the rocket is again untrimmed. 

The manner in which angle of attack is generated is shown in Figure 2, a vector 
diagram of the forces on the rocket immediately after launch. An^e BOC is the 
launch elevation angle. OC is the thrust minus drag vector, aligned with the 
pxincipal body axis. Vector CD represents the acceleration of gravity. Vector 
OD is the sum of OC (the thrust minus drag) and CD (gravity). The rocket is 
accelerated from rest in direction OD, but it is pointing in direction OC. Angle 
DOC represents the angle of attack. 

This angle of attack is gradually reduced to zero hy the stabilizing moment (Mo) 
produced by the rocket fins. Trimming does not take place immediately because 
the stabilizing tendency is resisted by the large moment of inertia about the 
rocket's principal axis. 

If the an^e of attack is assumed to be zero during Period I (i. e. , before the rocket 
trims out), large errors in the predicted trajectory may result. This is the case 
when standard three linear degree of freedom simulations are used. Further 
errors may result from assumptions about other initial and launch conditions. 

To reduce these errors, the computer program described in this paper uses three 
separate subsections to predict trajectories. A launch rail subsection is used 
until the rocket has left its launcher. The program then switches to a special 
'3-D' section which computes motions in two linear and one angular degrees of 
freedom. This permits accurate simulation of Period I flight when the angle of 
attach is finite. When the rocket trims out, the program switches to the standard, 
three linear degrees of freedom model. This model is used throughout Periods n 
andm. 



ANGLE OF ATTACK (a) 
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FIGURE 1; Angle-of- Attack/Time History 
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OC = Thrust-Drag 
■fCD = Gravity 
OD = Resultant 


Forces on Rocket 

immediately after Launch 
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The standard '3-D' model can be used accurately in Period III despite the un- 
trimmed condition because the rocket is not thrusting. This untrimmed condition 
results from the lack of restoring moment due to low atmospheric density. The 
aerodynamic force on the fins varies as the product of velocity squared and at- 
mospheric density. Thus, although the velocity is hi^ in Period m, the density 
is low enough that the restoring moment is insufficient to maintain a zero an^e 
of attack. 

In the absence of thrust, therefore, the untrimmed fli^t of Period m has no 
elKect on the trajectory. Should there be an upper stage to be fired in Period HI, 
the "HYBRID 3-D" program has a thrust vector control option to handle !his case. 
There is also an option to change the aerodynamic drag on the spent rocket during 
this period. 

^y using the launch rail subsection, the special '3-D' section, the standard '3-D' 
section, and the thrust vector control and aerodynamic drag options, the program 
can calculate the rocket's impact point with sufficient accuracy for range safety 
hazard > stimation and for the planning of payload recovery operations. 

B. BASIC ASSUMPTIONS 

The mathematical model used in the "HYBRID 3-D" program includes such effects 
as the variation of aerodynamic coefficients with Mach number, the change of 
engine thrust with time, and the change of rocket mass due to propellant burning 
and stage separation. But some simplifying assumptions have been necessary to 
speed calculation. These are detailed below. 

Please note that where equations are given, the notation is similar to that of 
FORTRAN, with the provision that all variables are of type 'IISAL'. 

1. Earth's Geometric Shg^pe 

The shape of the earth is a "pear" shaped spheroid with a sli^tly smaller 
northern hemisphere. The oblate spheroid is usually used as an approxi- 
mation to the earth's shape. However, for sounding rocket trajectories. 



where the range of the rocket is short in comparison to the circumference 
of the earth, a simplified model may be utilized: the "local spherical 
earth" in which we approximate the earth's shape by a sphere with a 
radius equal to the geocentric radius at the launch site. In the program, 
the local geocentric radius of the laimch site (RE) is set to 20899262. 
feet, a value appropriate for Wallops Island. 

2. Earth's Gravitational Field 

The gravitational potential of the earth is simplified to the inverse square 
law without harmonics, i. e. : 

GG a GM/R^+3 (B-1) 

where: GG - is the gravitational potential; 

GM - is defined as 1.4076576E16 feetxt*3/secc'nds*>t'2; 

R - is the distance from the geocenter to the vehicle. 

The components of the gravitational potential are: 

GX => -GG»X 
GY =< -GG*Y 
GZ — -GG*Z 

where X, Y, and Z are the geocentric coordinates of the vehicle. 

3. Atmosphere 

We have selected an atmospheric model which closely approximates the 
1962 U. S. Standard Atmosphere. (Reference 6) 

The subroutine ATMSPH is called with altitude as the argument. The 
atmospheric parameters pressure, density, temperature, viscosity, 
and speed of sound are returned. 

4. Earth's Rotation 


( 


During the short duration of Period I, the rotation of the earth may be 
neglected without appreciable error. In Periods n and m, the effect of 





the earth's rotation is accoimted for in a transformation from a Launch 
Inertial Coordinate System to an Earth Fixed Coordinate System. 

C. COORDINATE SYSTEMS 

In the trajectory calculation, five sets of Cartesian coordinates are used. They 
are discussed in detail below: 

1. Inertial Coordinate System (X, Y, Z) (Inertial at Launch): 

The origin of the Inertial System is the earth's center. Once the X and 

Y axes are determined at launch time, the Inertial cocrdinates are fibced 
in earth-centered space, and do not rotate with the earth. 

X - on the earth's equatorial plane, pointing to zero longitude at 
launch. 

Y - on the earth's equatorial plane, pointing to 90^ East longitude 

at launch. 

Z - perpendicular to the equatorial plane, pointing to the North Pole. 

2. Earth- Fixed Coordinate System (X , Y_, Z ): 

£ £ £ 

The origin of the Earth- Fixed coordinate system is the center of the 
earth. However, unlike the Inertial system, whose origin is also the 
earth's center, the X^ and Y^ axes of the Earth- Fixed system are the 
Greenwich and 90*^ longitudes respectively and rotate with the earth. 

At the moment of launch, the Earth- Fixed axes coincide with the Inertial 
axes. 

X^ - on the earth's equatorial plane, always pointing to the Greenwich 
longitude. 

o 

Y^ - on the earth's equatorial plane, always pointing to 90 East 
longitude. 

Z_ - perpendicular to the equatorial plane, pointing to the North Pole. 
£ 



3. Instantaneous-Topocentric Coordinate System (x,y,z> ! 

The origin of the Instantaneous-Topocentric coordinates is the projection 
point of the moving rocket on the earth's surface, the point at which the 
geocentric radius vector to the rocket intersects the earth's surface. 

X - on the local horizon plane tangent to the instantaneous projection 
of the rocket, directed along the local geocentric north. 

y - on the local horizon plane tangent to the instantaneous projection 
of the rocket, directed along the local geocentric east. 

z - perpendicular to the instantaneous local tangent plane, directed 
along the geocentric radius vector, and pointing toward the 
earth's center to complete the ri^t-hand system. 



The origin of the Launch coordinate system is the launch site. The 
axis is chosen to point in the direction of launch. 

- on the launch-tangent plane, pointing in the direction of launch. 

y_ - on the launch-tangent plane, pointing normal to the launch 

azimuth in the direction which with x_ and z forms a ri^t- 
band system. 

z^ - perpendicular to the launch-tangent plane, positive upward 
from the earth's center* 


ft 




5. Body Coordinate ^stem 

The origin of the Bo^ coordinate system is the center of mass of the 
rocket. 

X- - along the rocket principle axis, positive forward. 

y - normal to the plane in the direction which completes 

the lig^t-hand ^stm. 

z - perpendicular to the x_ axis and contained in the plane of 

ssrmmetry of the rockS, positive downward. ; 
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D. DYNAMIC EQUATIONS OF THE ROCKET 

••HYBRID 3-D" simulates the trajectory of the rocket with two models, one for 
Period I and another for the rest of the flight. The dynamic equations for eac'^ 
model are given here. A FORTRAN-like notation is used (all variables arf 
type "REAL"). 


1. The Equations for Period I, using two linear and one angular degrees 
of freedom: 


AXL =* ((THRUST-DRAG).t.CS - FNORM*SN)/MASS 

AZL « ((THRUST-DRAG)*SN + FNORM*CS)/MASS-GO 


THDD 

where 

AXL 

THRUST 

DRAG 


(-LSM*FNORM)/INERT 


acceleration in the XL direction, which is downrange at 
launch 

the thrust of the rocket motor. Found from table look-up, 
with time as the argument 

CDS*Q 


Q 

VL 

CDS 


Q 

DENS 


(DENS*VL*t«2)/2.0 
SQRT(VXL**2 4- VZL»*2) 

is given by sub- routine TABl which carters a look-iq) table 
with Mach number as the argument to find CD, the co- 
efficient of drag, and multiplies this by SAREA, the 
reference area 

is called the "dynamic pressure". 

Is given by the Standard Atmosphere sub-routine 'ATMSPH', 
entered with altitude as the argument. This sub-routine 


also gives 


TEMP 

temperature 

PRES 

atmospheric pressure 

Vise 

viscosity 

SOUNT 

the speed of sound. 


( 



VXL.VZL 

CS 

THETA 


ENORM 


are velocity components in the downrange and vertical 
directions. Found by integration of AXL and AZL. 

COS{THETA) 

is the inclination angle of the rocket, the an^e between the 
principal bo^ axis and the hciiizontal. It is originally set 
to the launch elevation angle, and is modified by double 
integration of THDD. THETA double dot. 

CNA.«<SAREA*Q»ALPH 


CNA 

SAREA 


is the slope of the coefficient of the normal force acting 
on the center of pressure. 

is the reference area. 


Q • is defined above. 


ALFH 


GAMA 


MASS 

AZL 


THETA > GAMA, and is called the **an^e of attack*'. The 
program switches from the Period I model to the model 
used for the rest of the flis^t when ALPH becomes zero 
(found by interpolation as it crosses from plus to minus). 
The duration of Period I is about 0. 8 second for a 
NIKE-CAJUN. and about 5 seconds for a SCOUT. 

ATAN (VZL/VXL) 

GAMA is the flig^ path angle which defines the dii . 
in which the rocket is moving 

is the wei^t of the rocket divided fay the force of gx..ity. 

is acceleration in the Z direction, which is vertical at 
the launch site. “ 


THDD - is THETA double dot. the acceleration in the rocket's 
inclination angle to the axis. 

LSM a LCP-LCG. 

LCP • is the distance from the reference position (usually the 
nose of the rocket) to the center of pressure. 

LCG - is the distance from the reference position to the rocket's 
center of gravity. 

INERT - is the rocket's pitch moment of inertia. 


The equations above have been simplified Ity the elimination of efiects such 
as pitch damping and Jet damping which were found to have no significant 



effect on the solution, since Period I is of such short duration. Some 
variables are made constants, with the same justiff cation. Some of 
these are: GO; CNA; LSM; and INERT. 


2. The Equations for Periods II and IH: 

The angular motion of the rocket is neglected; only the linear motion 
is considered. The equatior'i of motion are written in the Inertial 
Coordinate System. 


XDDOT - 
YDDOT - 
2DDOT 
XDDOT - 

GX 

THX 


UX 

MASS 

YDDOT, 

ZDDOT 


GX + (THX - UX)/MASS 
GY + (THY - UY)/MA.SS 
GZ + (THZ - UZ)/MASS 

is X double dot; the acceleration is the direction of the 
Inertial X-axis (please i*efer to Section C for a definition 
of the various coordinate systems used). 

is the X component of the acceleration of gravity, aa 
felt at the rocket's center of gravity. 

is the rocket's component of thrust in the X-direction. 

This component of thrust would nozmally be taken in the 
direction of rocket motion, except that it may be modified 
by the thrust control option. The thrust control option 
may be specified to hold the thrust direction of the rocket 
constant after r giveu time of flight or rocket altitude. 

This option is useful when the rocket attains an altitude 
where the atmospheric density is so low that the rocket's 
gyroscf^ic stability tends to keep i!: pointing in a constant 
direction and there is to be an uppe-stage firing (Period UI). 

is the component of drag in the inertial, X direction. 

is the present weight of the rocket, derived from weight 
tables, divided by the sea-level acceleration due to gravity. 

are the components of acceleration in the inertial Y and Z 
directions. Their equations are symmetrical with the one 
for XDDOT. 
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3. The Equations for the Launch Rail: 

The user may specify the type of latmcher to be used for the simulation. 
The simulated "ocket will then be accelerated from rest along the launch 
rail and will have an initial velocity at release. This initial velocity 
will decrease the duration of Period I. The equation of motion is: 


LACC 
LACC 
THRUST - 
WEIGHT - 


(THRUST - WEIGHT*SEL - DRAG)/MASS 

is acceleration in the direction of the launch rail 

is found by table look-up, with time as the argument 

is the rocket weight, found by table look-up; time is the 
argument. 


SEL - is the sine of the elevation angle of the laimcher. 

DRAG - is the aerodynamic dr^ of the rocket. This term is quite 

small, and could be neglected. Tae friction drag of the 
launcher is not included, although the equation might be 
reformulated to add it to the aero(fynamic drag. In the 
practical case of launch elevation angles of perhaps 80^, 
the normal force on the launch rail (which is multiplied 
by the coefficient of friction to give the friction drag) is 
very low, resulting in negligible friction drag. 


MASS a WEIGHT/GO 

GO - is the acceleration due to gravity at sea level. 


The initial conditions are entered and then this equation is integrated 
repeatedly to give the velocity and distance traveled on the launch rail. 
When this distance equals the length of the laimcher, the program switches 
to solution of the equations for Period I. 
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SUMMARY 


A computer program, "HYBRID 3-D", has been described which simplifies tra- 
jectory simulation. 

A comparison of a "6-D" trajectory, a "standard" 3-d trajectory, and "HYBRID 
3-D" is shown in Figure 3. "HYBRID 3-D" achieves good accuracy through the 
special treatment of the untrimmed period near laimch. 

"HYBRID 3-D" also contains a simulation of the dynamics of the launch rail to 
enhance accuracy in simulating the trajectory of slow-accelerating rockets. 
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INTRODUCTION 


This trajectory program is written in Fortran IV and was tested 
on the Honeywell 625 computer. Input data is transferred from 
punched cards to disc (file code 05). The program reads this 
disc, prints the data, rewinds the disc, and rereads the disc 
as data is needed for calculations . , The math model includes 
aerodynamic drag, rotating earth, launch rail simulation, rigid 
body dynamics at lift-off, and standard 3D point-mass equations 
After leaving the launcher, the rocket is restricted to motion 
in two dimensions until the angle of attack goes to zero; then 
the motion becomes three dimensional with thrust and velocity 
vectors aligned. The equations of motion are Integrated by the 
Runge-Kutta method, utilizing variable step size for better 
efficiency. Processing time is approximately 0.7 minutes for 
each 100 seconds of trajectory time. Output data may be 
printed in English or metric units, and written on tape (file 
code 11) . 



2 . 


PROGRAM VARIABLE DEFINITIONS 




( 


( 


A 

ACX,ACY,ACZ 

ACCI 

ALIM 

ALPH 

ALT 

ALTC 

AREA 

ARRAY 

ARREY 

AXL 

AZ 

AZL 

AZLCH,ELLCH 

AZRAD,ELRAD 

CAPT 

CD 

CDD 

CNA 

CONV 

CPHI 


semi-major axis of IIP ellipse 

components of acceleration with respect to 
topocentrlc system 

inertial acceleration 

altitude where pressure and density are set 
equal to zero (400,000 ft.) 

angle of attack (o) 

altitude 

input altitude at which thrust vector remains 
in constant direction 

rocket nozzle exit area (Aj^) 

matrix for storing all output data (in English 
units) 

matrix for storing all output data in metric 
units 

acceleration along launch azimuth in launch 
coordinate system 

flight azimuth angle; or launch azimuth (A) 

acceleration in vertical direction in launch 
coordinate system 

look azimuth and look elevation angle with 
respect to launch site 

look azimuth and look elevation angle with 
respect to radar 

time from launch to Impact at IIP 

drag coefficient (C^) 

input table of drag coefficients (Cj.); 
(Includes drag area, if SAREA » 1.0; 

normal force coefficient 

if/l80® 

cos (♦) 







DELTA 


DENS 

DIPP 

DIST 

DRAG 

DT 

EDLCH,NDLCH, 

ZDLCH 

EDRAD,NDRAD, 

ZDRAD 

EL 

ELCH,NLCH, 

ZLCH 

ELPREV 

EPBIG 

EPS 

EPTINY 

ERAD,NRAD, 

ZRAD 

PNORM 

PRAC 

GAMA 

GO 

GDLAT 


■ range angle (6) 

» atmospheric density 
» thrust minus drag 
“ distance traveled on launch rail 
» drag force (D) 

* Integration time Interval 


® east, north, vertical velocity components 
with respect to launch site 


= east, north, vertical velocity components 
with respect to radar 

= flight elevation angle or launch elevation 
angle (r) 


= east, north, vertical coordinates with respect 
to launch site 

= previous value of velocity elevation angle 

= maximum relative error allowed in Runge-Kutta 
method (1.0 x 10“5) 

= geodetic latitude minus geocentric latitude 

= minimum relative error allowed in Runge-Kutta 
method (5.0 x 10"7) 

= east, north, vertical coordinates of rocket 

with respect to radar 

= normal force on rocket 

= Interpolation fraction for Mach number 

■ velocity elevation angle (r) 

■ gravity acceleration at Earth surface 
( 32.174 ft/sec2) 

« geodetic latitude 





( 


) 

% 


( 


GM 

GX, GY, GZ 
HIGH 

lEND 

IIPLAT,IIPLON 
IIPR,IIPTIM = 


1 n O 

gravitation constant (1.^076576 x lO"*" ff^/sec ) 

components of gravity acceleration 

maximum integration interval ■ one half the 
print time interval 

control constant (set « 1 at Impact) 

latitude (GD), longitude, ground range, and 
impact time for instantaneous impact point 


INERT 
INEXT 
I OPT 

IPRINT 

K 

KCON 

KDEL 

KPAG 

KSPENT 

KSTOP 

LACC 

LAMDA 

LAMDAO 

LAMDAl 

LAT 

LENGTH 

LINE 


® moment of inertia of rocket (I) 

= control constant (set = 1 at each phase time) 

* thrust reference option (=1 for sea level, 

=2 for vacuum) 

= control constant (set = 1 when printout is 
desired) 

= matrix for storing K values for Runge-Kutta 
method 

= 1 for main rocket calculations; 

2 for spent stage calculations 

= Increment for subscript of phase control 
constant 

= output page number 

= present spent stage number 

= 0 for no stop at apogee 

1 for stopping calculations at apogee 

= acceleration on launch rail 

= angle from X axis in X, Y plane (X) 

® launch longitude (Xq) 

» radar longitude (x^) 

= geocentric latitude ((>) 

= length of launch rail 

» number of printout line 


i 






LONG 


longitude (X) 


LOW 

m 

smallest integration time Interval 
(1.0 X 10-4 sec) 

LSM 

m 

distance from center of gravity to center 
of pressure (n_,„) 

LVEL 

m 

velocity on launch rail 

MACH 

s 

Mach number (M) 

MASS 

s 

total rocket mass - 

MCH 

s 

input table of Mach numbers corresponding to 
table of drag coefficients 

MESS 

s 

array for storing printout message 

NAP 

s 

control constant (= 1 at apogee) 

NCD 

8S 

number of drag coefficients in table 


NDEL = Increment for subscript of phase time 

NGEO = 0 for geocentric elevation angle and azimuth 

(input) 

1 for geodetic elevation angle and azimuth 
(input) 

NID « control constants for phase changes 

= 1, read phase message only 
« 2, read phase message and drag table 
» 3, read phase message and thrust, weight, 
drag tables 

= 4, nothing is read in 

NIIP « input constant for special output tape 

■0, tape time interval same as print time 
interval 

■ 1, tape time interval equal 0.1 sec, and 

PTI must be an integer multiple of 0.1 
sec 

■ 2, no tape output 

NMT ■ nximber of rocket motors 


NPAGE “ 1 for output page A 

2 for output page B 

3 for output page C 

4 for output page D 

5 for output page E 

6 for output page F 



i 
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NPH » 
KPHS - 
NSEP - 
NSKIP - 

NSPENT = 
NSTOP = 

NSYS » 


NTH » 
NWT = 
OMEGA = 

PO = 

PHAS = 
PHI « 
PlilO. = 
PHIl « 
PHT « 

PI « 

PI2 » 
PRES « 
PTI 

PTIME • 


present phase number 

number of phase times 

number of motor separations 

input control constant 
« Of ALPHA routine performed 
• 1 , ALPHA routine skipped 

number of spent stages left to be calculated 

0 for TSTOP ® 0.0 (no termination on-time) 

1 for TSTOP / 0 (trajectory terminates on time) 

output control constant 
1, for English units 
2f for metric units 
3f both English emd metric 

nuniber of thrust values in table 

number of propelleuit weights in table 

earth rotation speed ( 7,29211 * 10“^ rad/sec) 

(«) 

atmospheric pressure at surface of earth 
(2115.666 lbs/ft2) 

output message in column 2 to 7 

angle from X, Y plane to R vector (^) 

launcher geodetic latitude 

radar site geodetic latitude 

phase times (input) 

3.14159265 (ir) 

2ir 

atmospheric pressure (p) 
print time Interval 
print time 


Q - 

R - 
RA ■ 
RANGE « 
RB » 

RE “ 
RLCH B 
RRAD s 
S e 

SAREA = 

SA « 
SMAT a 
SOUND = 
SPHI » 
SPS a 
STPTI = 
TO a 
T1 a 

TEMP a 
TFRAC a 
THCON a 


THDD a 
THETA - 
THREE ■ 


dynamls pressure 

distance from earth center to rocket 

equatorial radius (20925741 ft) 

ground range from launch site to rocket 

polar radius (20855591 ft) 

earth radius (20899262 ft) 

slemt range from launch site to rocket 

slant range from radar to rocket 

matrix for storing X, Y, Z, X, Y, Z, X, Y, Z 
in Runge-Kutta routine 

drag area (* i.o when area is included in 
drag coefficient) 

matrix for storing all output data 
transformation matrix, topocentric to inertial 
speed of sound 
sin (4) 

array for storing initial spent stage variables 

stored value of the print time interval 

launch time » 0.0 sec (t^j) 

time to begin constemt thrust dlxrection in 
inertial system 

atmospheric temperature 

time interpolation fraction 

1 for no thrust control 

2 for constant thrust direction beginning at 
time Tl, or at altitude ALTC 

angular acceleration in launch coordinate system 

thrust elevation angle 

thrust at sea level or in vacuum 



I 




THRUST « 

rocket thrust (TH) 


• 

THS « 

thrust table input 



TIM - 

time table corresponding to thrust table 

’ 


TIMl, TIM2 * 

times used to control writing output tape 

, 


TIME ■ 

time from launch (t) 

i 


TMAT * 

transformation matrix, inertial to topocentric 

1 

j 

i 

I 


TSEP « 

table of separation times rc x single stage 
rocket, set TSEP » 0.0) 

1 


TSTOP « 

time to stop calculations 

I 

1 

1 


TWT * 

time table corresponding to propellant 
weight table 



VEL s 

total speed of rocket (inertial) 

! 


VELLCH « 

rocket speed with respect to launch site 

i 


VELRAD B 

rocket speed with respect to radar 


( 

Vise ■ 

atmospheric viscosity 

; 


VT « 

total velocity (topocentric) 



VX, VY, VZ » 

topocentric velocity components 



WEIGHT « 

total rocket weight 



WPL * 

payload weight 



WPP •* 

present propellant weight 



WPP2 ■ 

total propellant weight of all unfired motors 



WPR ■ 

input table for propellant weight 



WRM > 

present rocket motor weight 

4 

I 


WTM * 

table of rocket motor weights 



WTP « 

initial propellant weight in each rocket 

i 


WTS - 

spent stage weight 



X, Y, Z - 

inertial coordinates 


I 

j. 


5 



XDOT, YDOT, 
ZDOT » 


inertial velocity components (X, Y, Z) 


XDDOT, YDDOT, 
ZDDOT - 

XIP, YIP, ZIP - 
XL, YL, ZL = 

XT, YT, ZT - 


inertial acceleration components (X, Y, Z) 

inertial position of the instantaneous impact 
point 

launch system coordinates where XL is along 
launch azimuth 

topocentric coordinates 



3. RUNGE-KUTTA INTEGRATION ROUTINE 

For a discussion of the Runge-Kutta formula of order 4, see 
Reference 7. 

Here, only the equations in X are listed, since the Y euid Z 
equations correspond exactly. N 

(NOTE • the actual subroutine uses matrix S to store all the 
inertial variables, X, Y, Z, X, Y, Z, X, Y, Z) 

• ** 

Step 1: Input Xq, Xq, Xq,* time ** T^ 


"ll- 

' (DT) Xq 

K41 - 

(DT) 

*0 

tire 

= Tj^ = T 

0 + DT/2 



='i = 

Xq +(^) 

^0 



■ 

Xq 

^0 



*1- 

X (Tj^, X 

1' ^1' ^1 



*12 ■ 

‘ (DT) 

K42 = 

(DT) 

*1 


' m 









DT 


Step 3: time * * Tq + j— 


*2 " *0 * <“> *1 


• • TYV 

*2 * =‘o ♦ 'P \ 


(NOTE! T, • Tj) 


Xj . X (Tj, Xj, Yj, Zj, Xj, Yj, Zj) 
Kj 3 . (DT) ij K^3 . (DT) ij 


Step 4: time = * Tq + DT 

X 3 = Xq + (DT) X 2 
X 3 = Xq + (DT) X 2 

*3 = ^3' ^3' ^3» *3» ^3' V 

» (DT) X 3 = (DT) X 3 

Solution at time » T- * T« + DT 

4 0 


X» Xq + 

(1/6) 

(K33 + 

2Ki2 + 

2X^3* 

*^14 

i » io + 

(1/6) 

(Ki, + 

2K42 + 

2K43* 

^44 

X = X(T^, 

X, Y, 

Z, X, 

Y, Z) 




Error Analysis : 


The Runge-Kutta method of order 4 will result in an error of 
order 5, which la reduced by using "extrapolation to the 
limit" . 


First, find a solution using DT » H; then find a corresponding 
solution by using DT s (H/2) twice. Finally, combine the two 
solutions in such a way that most of the order 5 error is 
eliminated. 

Xgjj * exact solution 

X(l) = solution using DT = H 
X(2) * solution using DT « H/2 twice 
X(l) = Xgj^ + + order 6 

% X(2) = Xgjj + B(H/2)^ + C(H/2)® + order 6 

The factors A, B, and C are composed of derivatives of the 
acceleration function. The X(2) solution has two order 5 
errors, because twice as many steps are required when 
DT « (H/2) . 

If H is small. A, B, and C will be approximately equal, 
resulting in: 

X(l) - Xgjj : AH^ 

AH^ 

X(2) - Xgjj : ^ 

Xgjj : X(2) + (1/15) [X(2) - X(1)1 

The right hand side is an improved approximation for the exact 

• • • 

solution. Similar equations are used for X, Z, X, Y, Z. 


Optimum DT 


In the Runge-Kutta formula, the theoretical error term is of 
order 5, and the smaller the step size, DT, the smaller ihis 
error becomes. However, when DT is made smaller, the nximber 
of integration steps increases, resulting in more computer 
round-off error. 

To find an "optimum DT", the program proceeds as follows: 

• •• • 

1. Find the velocity components X(l), Y(l), Z(l), using 
DT = H. 

2. Find X(2), Y(2), Z (2) , using DT = H/2 twice. 

3. Calculate the relative errors in velocity: 

E = .1^,(1); - I , etc. 

X(2) 

4. If E < EPTINY, then DT is too small, and computer time 
would be wasted. Thus, DT is increased for the next 
integration interval. 

5. If E > EPBIG, then DT is too large, and the velocity 
error is unacceptable. Thus, the calculation is repeated 
using a smaller DT. 

The choice of EPTINY and EPBIG is quite arbitrary. 

Good results have been obtained with EPTINY = lO”^ 
and EPBIG = 10"^. 



4. general plow chart 
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TIME TO STOP CALCULATIOIIS 
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5. SUBROUTINE DESCRIPTION 


ACCELO Subroutine ACCELO computes the relative linear 

acceleration and angular acceleration, assximing 
the rocket flies in the launch plane (2 dimensions) 
with an angle of attack. 

ACCEL 1 Subroutine ACCEll computes the tl-iree components of 

inertial acceleration which determine the main 
trajectory. 

ACCEL2 Subroutine ACCEL2 computes the relative linear 

acceleration (one dimension) for motion along the 
launch rail, assuming there is no frictional force. 

ACCEL3 Subroutine ACCEL3 computes the acceleration using 

elevation and azimuth rather than velocity com- 
ponents . 

ALFA Subroutine ALFA determines the trajectory while the 

rocket flies with an cuigle of attack during the launch 
phase. Linear motion is constrained to the launch 
plane until the angle of attack goes to zero (or until 
the phase is nearly completed) . 

^ ATMSPH Subroutine ATMSPH computes temperature, pressure, 

; viscosity, density, and speed of sound using an 

I 

eighth degree polynomial approximation for the 
1962 Standard Atmosphere data. 

r 


I 
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CNALPA Subroutine CNALFA interpolates to find the normal 

force coefficient and pitch damping coefficient 
which are used in the angle of attack equations. 

DATAl — - Subroutine DATAl computes most of the output data, 

DIRECT Subroutine DIRECT determines the direction of the 

thrust vector; the thrust is either lined up with 
velocity, or the thrust maintains a constant di- 
rection with respect to the inertial system. 
Normally, the first option is used; however, for 
spin stabilized rockets, the second option can be 
used at high altitudes. (NOTE: This subroutine is 
not called during the ALFA routine.) 

IIP — Subroutine IIP computes the instantaneous impact points 
assuming zero thrust emd vacuum trajectory. 

INPUT — - Subroutine INPUT reads data from input tape 05 and 
prints it. 

LAUNCH — - Subroutine LAUNCH simulates a frictionless launch 

rail tut includes aerodynamic drag in the equations 
of motion. 

MAT Subroutine MAT calculates the elements of the trans- 


formation matrix for a coordinate system rotation 
The inverse matrix is also calculated. 




MATLCH — > Subroutine MATLCH calculates the transformation 

from the inertial system to the launch site system. 

MATRAD Subroutine MATRAD calculates the transformation 

from the inertial system to the radar site system. 

MATROC Subroutine MATROC calculates the transformation 
from the inertial system to the instantaneous 
topocentric system. 

OUT Subroutine OUT interpolates all output data, stores 

it, and prints it out. 

RUNGE Subroutine RUNGE integrates the equations of motion, 

using the fourth order Runge-Kutta method with 
variable step-size. Refer to Part 2, Section 3. 

SETSA Subroutine SETSA transfers output data to a temporary 

storage array. 

TABLE — Subroutine TABLE reads in thrust, weight, emd drag 
tables and phase messages. 

TABl Subroutine TABl interpolates to find the drag coefficient. 

TAB2 - — Subroutine TAB2 interpolates to find the thrust. 

TAB3 — Subroutine TAB3 interpolates to find the total weight. 

TPOUT Subroutine TPOUT writes a special output tape with 

time interval controlled by NIIP. (See input 
description. ) 


I 


TRANSl — Subroutine TRAKSl trwsforms from the inertial 


system to one of the rotating systems. 

TRANS2 Subroutine TRANS2 tremsforms from a rotating sys 

tern to the inertial system. 
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HYBRID 3-D INPUT DESCRIPTION 

Input cards must be in FORTRAN NAMELIST format or in the form 
of a message (see sample); all cards must be punched only In 
columns 2 through 72. Each input value must be followed by a 
comma, except the last item in a list, which is followed by a 
$. 

» 

1. The following lists must always be input: NAMQ, NAMR, 

NAMS, NAMT, NAMU, NAMV. Within a list, if a particular 
variable is not required, it should be omitted. 

2. The lists NAM2, NAM3, NAM^, NAM5 may be input depending 
upon the trajectory and type of rocket. Thrust, weight, 
and drag table input is controlled by the phase time 
array, PHT, and the array NID. 

3. All thrust, weight, and drag tables have a limit of 

50 values and must have at least two values. The input 
arrays PHT and NID have a limit of IM values; the 
arrays WTM and WTP have a limit of 6 values. 

4. For a thrusting phase, the thrust, weight, and drag tables 
are preceded by one card containing a "phase message"; 
this may be any phrase to identify the phase; e.g. STAGE 

1 THRUSTING. A phase message card must also precede the 
drag table for a coasting phase. The phase messages may 
be punched in coliunn 2 through 72, but only the characters 
in column 2 through 7 will appear on the output (in column 

2 through 7) (See sample input and output.) 


i 
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5. To simulate the rocket launcher, set LENGTH equal to 
the length of the launch rail. This option restricts 


the motion to one dimension only. 

6. To simulate angle of attack near lift-off, set NSKIP ■ 0. 
Motion will be In two dimensions until the angle of 
attack goes to zero; then, the program switches to 
three dimensions, with thrust and velocity vectors 
aligned. (Angle of attack can be calculated only during 
first stage thrusting.) 

7. If NSKIP ■ 1, then the angle of attack routine Is 
skipped, and the motion will be In three dimensions with 
thrust and velocity aligned. 

8. The trajectory terminates (1) at altitude equal zero; 

(2) at apogee; (3) or at a given time. See the control 
constants NSTOP, TSTOP, KSTOP. If the trajectory Is not 
stopped at apogee or at a given time, then It terminates 
when altitude equals zero. 

9. For tape output, see control constant NIIP. In the deck for 
execution, include a tape control card with file code 11. 



I 


INPUT DESCRIPTION 


FIRST CARD; title phrase in column 2 to 72 


NAMQ List : 


$NAMQ in column 2 through 6 followed by; 


VT 

EL 

AZ 


initial speed (ft/s-c) 

initial flight elevation angle or launch elevation 
angle ((’eg.) 

initial flight azimuth angle or launch azimuth angle 
(deg. ) 


NGEO 

PHI 

LAMDA 

TIME 

ALT 

LENGTH 


0 if EL, AZ are geocentric angles 

1 if EL, AZ are geodetic angles 

initial geodetic latitude (deg.) 

initial longitude (deg.) 

initial time (not less than Q.O), (sec) 

(launch time should equal zero) 

initial altitude (ft.), (NOTE: launch altitude should 
equal 0.0 only) 

length of launch rail (ft.) 


NAMR List: 


$NAMR in column 2 through 6 followed by; 


PTI 

NS PENT 

NSTOP 

TSTOP 

NSKIP 

KSTOP 


print time Interval (not less than .01), (sec) 

nviraber of spent stages to be run 

1 if TSTOP not equal 0.0 (otherwise, omit) 

time at which calculations are terminated (sec.), 
(if not desired, omit) 

0 alpha routine calculations (beginning at launch) 

1 alpha routine skipped 

1 to stop calculations at apogee (otherwise, omit) 


THCON ■ 1 for no thrust control 

2 for constant thrust direction beginning at time ■ 
T1 or at altitude ALTC (the constant direction is 
determined by the program) 


T1 

ALTC 

NSYS 

NIIP 


NPAQE 


time to begin constant thrust direction (sec.), 

(if not used, omit) 

altitude where thrust co;.trol begins (ft.), (if not 
used, omit) 

1 for English output 

2 for metric output 

3 for both English and metric output 

0, special tape time interval same as print time Interval 

1, tape time Interval equal 0.1 sec., and PTI equal an 
integer multiple of 0.1 sec. 

2, no tape output 

1 for output page A 

2 for output page B 

3 for output page C 

4 for output page D 

5 for output page E 

6 for output page P 

Ce.g., if NPAGE • 2, 4, 5j pages B, D, E will pe 
printed out] 


NAMS List : 

$NAMS in column 2 throug^i 6 followed by: 

PHIO ■» launcher geodetic latitude (deg.) 

LAMDAO « launcher longitude (deg.) 

PHIl ■ radar site geodetic latitude (deg.) (if not used, omit) 

LAMDAl ■ radar site longitude (deg.) (if not used, omit) 


NAMT List ; 

$NAMT in column 2 through 6 followed by; 

NPHS/PHT- phase time array (sec.) (e.g. 1st stage ignition time, 
1st stage burnout time, 2nd stage ignition, 2nd stage 
burnout, etc.) 



NTAB/NID 


array of control constants for phase changes 


1, to read in phase message only 

2, to read phase message and drag table 

3, to read in phase message and thrust, weight, 
and drag tables 

4, nothing is read in 

[NOTE: There must be a value of NID for each value 

of PHT, plus one additional value corresponding 
to apogee.] 

NAMU List : 

$NAMU in column 2 through 6 followed by : 

NMT/WTM = list of inert rocket motor weights (lbs.) (Includes 
miscellaneous weight) 

WPL * payload weiglit (lbs.) 

NSEP/TSEP = table of motor separation times (sec.) 

[NOTE: For single stage rocket, omit] 

[Each value of TSEP must appear in the PHT array.] 

WTP = list of propellant weights (one weight for each stage), 

(lbs.) 


NAMV List : 

$NAMV in column 2 through 6 followed by: 

LSM = distance from center of gravity to center of 

pressure ( ft . ) 

2 

INERT “ moment of inertia (slug ft. ) 

CNA = normal force coefficient (rad~^) 


NAM2 List : (thrust table, input if NID = 3) 

.* $NAM2 in column 2 through 6 followed by : 

t JL AREA = rocket nozzle exit area (ft. ) 

lOPT = 1 for sea level thrust table, 

2 for vacuum thrust table ■ 

TIM » time table corresponding to thrust table (sec.) 

NTH/THS = thrust table (lbs.) ( 

I 



i 


I 

i 
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NAM3 List ; (propellant weight table, input if NID » 3) 

$NAM3 in column 2 through 6 followed by: 

TWT = time table corresponding to weight table (sec.) 

NWT/WPR = propellant weight table (lbs.) 

NAM4 List : (drag data, input if NID = 2 or 3) 

$NAM4 in column 2 through 6 followed by: 

2 

SAREA = drag area (ft. ) 

(=1.0 when the area is Included in the drag coefficient) 

MCH = table of Mach numbers corresponding to drag coefficient 

table 

NCD/CDD = table of drag coefficients 

NAM5 List : (input if NSPENT > 0) 

$NAM5 in column 2 through 6 followed by: 

PTI = print time interval for spent stage printout (sec.) 

WTS = spent stage weight (lbs.) 

NTAB/NID = array of phase control constants for spent stage 
trajectory : 

1, to read phase message only 

2, to read phase message plus drag table 
4, nothing is read in 

NOTE: there should be one or two values for NID — 

one value for the start of the spent stage 
trajectory and a second value for apogee. If 
the trajectory begins after apogee, then only 
one value is needed. 

Last card of input data : 

Z 99999 in column 7 through 12 


WM lii 


f 






7. HYBRID 3-d output DESCRIPTION 

Printed output can be in English and/or metric units by setting 
the input constant NSYS. In addition, there are six different 
page options which can be selected using the array NPAGE. Tape 
output is in English units only. Printout peculiarities are 
listed below: 

1. If the launch rail simulation is used (when LENGTH > 0), 
output will be suppressed while the rocket moves along 
the rail. 

2. At phase changes, there is double printout plus a phase 
message from the input data. 

3. Built-in messages include: 

TSTOP , printed when the trajectory is terminated at a 
given time. 

APOGEE , printed out near the actual apogee. 

ALT=0 , printed at the end of the trajectory when the 
altitude equals zero. 

IIP output is suppressed for spent stage trajectories. 

5. Above the atmosphere limit (400,000 ft.), the MACH 
number will be blank. 


Notes on Metric Units : 

nt = newton lbs. = pounds 1 ft. » 0.3048 m 

kg = kilogram ft. « feet 1 N.M. = 1.852 km 

m « meter N.M. « nautical mile 1 lb. » 4.4482 nt 

1 slug = 14.594 kg 
1 Ib./ft^ » 47.880 nt/m^ 




PAGE A: 


» 


( 


TIME 
ALPHA 
TH EL 
FL EL 
FL AZ 
ALT 
RANGE 
LAT GD 
LONG 

PAGE B ; 

TIME 
THRUST 
WEIGHT 
DRAG 
MACH 
DYN PR 
REL ACC 
REL VEL 
MASS 


time from launch (sec.) 
angle of attack (deg.) 
thrust elevation angle (deg.) 
velocity elevation angle (deg.) 
velocity azimuth angle (deg.) 
altitude (ft.) 
surface range (N.M.) 
geodetic latitude (deg.) 
longitude (deg.) 


time from lavinch (sec.) 
rocket thrust (lbs.) 
rocket weight (lbs.) 
aerodynamic drag (lbs.) 

Mach number (suppressed when altitude > 400,000 ft.) 

dynamic pressure (lb. /ft. 2) 

acceleration relative to earth (ft./sec,^) 

velocity relative to earth (ft. /sec.) 

mass of rocket (slugs) 


. PAGE C; Launch Site Coordinate System 


TIME 

XL,YL,ZL 

VXL,VYL, 

VZL 

RXY 

6AML 

VEL-L 

RLDOT 


time from launch (sec.) 
position coordinates (ft.) 

velocity components (ft. /sec.) 
tangent plane range (ft.) 

velocity elevation angle with respect to tangent 
plane (deg.) 

velocity with respect to launch site (ft. /sec.) 
rate of change of slant range 


( 


PAGE D: Radar Site Coordinate System 


TIME 

SL RANGE 
LOOK AZ 
LOOK EL 
VEL 
EAST 

NORTH 

VERT 

VEL-E 

VEL-N 

VEL-V 


time from launch (sec.) 

slant range from radar (N.M.) 

look azimuth (deg.) 

look elevation angle (deg.) 

velocity relative to radar (ft. /sec.) 

distance to the East (+) or West (-) in the tangent 
plane (N.M.) 

distance to the North (+) or South (-) in the tangent 
plane (N.M.) 

distance (+) above the tangent plane (N.M.) 

Ecist (+) or West (-) component of velocity (ft. /sec.) 
North (+) or South (-) component of velocity (ft. /sec.) 
vertical component of velocity (ft. /sec.) 


PAGE E: Instantaneous Impact Points (Vacuum Trajectory) 

[NOTE: Page E is suppressed for spent stages] 


TIME 
LAT GD 
LONG 
RANGE 
IP TIME 


time from launch (sec.) 

geodetic latitude of IIP (deg.) 

longitude of IIP (deg.) 

range of IIP from launch site (N.M.) 

time from laionch to instantaneous impact (sec.) 


PAGE F : 

TIME 

R 

VEL 

ACCEL 


Inertial Coordinate System 

B time from launch (sec.) 

» distance from earth center (ft.) 

= inertial velocity (ft. /sec.) 

= inertial acceleration (ft./sec.^) 



SPECIAL TAPE OUTPUT 


This tape is high density (800 BPI) and contains 71 words per record, as listed 

below (those not listed are equal to zero): 

Word Number 

Contents (Enalish Units! 

1 

TIME 

2 

WEIGHT 

3 

THRUST 

4 

VEL (inertial) 

5 

R (inertial) 

6 

ALT 

% 

7 

RANGE 

8 

MACH 

9 

Q 

10 

LAT GD 

12 

LONG 

14 

ALPHA 

18 

FLEL 

20 

FLAZ 

23 

THEL 

37 

REL VEL 

55,56,57 

XL, YL, ZL 

58,59,60 

VXL, VYL, VZL 

61 

RXY 

62 

GAML 

63 

VEL-L 

64 

RLDOT 

65 

IP TIME (HP) 

66 

LAT GD (HP) 

67 

LONG (HP) 

68 

RANGE (HP) 


"i 
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Deck Set-Up with Object Program on Tape 


$ IDENT 153200 , HYBRID 3D 
$ EXECUTE 

$ LIMITS 100, 30K, 0, lOK 

$ TAPE R*, X1D„ 3969 

$ DATA 05 

$ INCODE IBMP 


(data cards) 

$ ENDJOB 
***EOP 


Deck Set-Up with Object Program on Cards 


$ IDENT 153200 , HYBRID 3D 
$ OPTION FORTRAN 

(object deck) 

$ EXECUTE 

$ LIMITS 100, 30K, 0, lOK 
$ DATA 05 

$ INCODE IBMP 


(data cards) 

$ ENDJOB 
»*»EOF 



9. 


SAMPLE INPUT AND OUTPUT 


I 


This Is a Nike-Cajun trajectory simulation using a 2-lnch 
launch rail and initial velocity equal to zero. Pour phases 
are shown: stage 1 thrusting, stage 2 coasting, stage 2 

thrusting, and stage 2 coasting. Input tables must be in the 
order specified by the NTAB/NID array: 

3 for stage 1 thrust; weight, drag input (NAM2, NAM3, NAM4) 

2 for stage 2 drag input for coasting (NAM4) 

3 for stage 2 thrust; weight, drag input (NAM2, NAM3, NAM4) 
ite 2 for stage 2 drag input for coasting (NAM4) 

4 for no input at apogee 


( 


4 

i 

/ ’ 



f 


1 


i 




The beginning of each phase is given by the NPHS/PHT array: 

-0.017 sec for stage 1 ignition 
3.523 for stage 1 burnout (and separation) 

17.0 for stage 2 ignition 

22.0 for stage 2 burnout 

Note that the times used for the thrust and weight tables are 
relative times which are added to the phase times by the pro- 
gram. If the first stage thrust table begins with a value less 
than the lift-off weight, then the ignition time must be adjusted i 
so that lift-off will occur at zero time. • 

s 

i 

1 

I 

5 
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1 
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NIKP-CAJUN TRAJPCTORY SIMULATION _ 

$mA MO V T = 0 • 0., a0« 0« A 2 _= J_20 . n_» NGFO = !♦ 


PHI s 37.flAB0» LAMO A = - 75,4736* TIME = n.O» ALT. =0.0* 
length =_.1667$ ; 


$MAMR 

PTI =0,5, NSPENT = 0, NSTOP = l* -fsTOP = 25. Q 

, ^ NSKIP=0* 



KSTOP * 0* THCCN = 1, Tl=Q.O* AXTC=0.6* 




nSYSsI , MIIP=0, MOAGEa 1,2$' 



$mAMS 

PHIO = +37.8480, LA'^OAO = -75.4736, 




DHIl = +‘^7.8412, LAMOAl = -75.4855$ 



SWAMT 

MPHS/oHT = -.017, 3.523, 17,0, 22.0, 




N'TAB/NIO = 3 , 2, 3, 2, 4$ 



SmAMU 

NMT/WTM = 564.0. 07.4, WPL = 50.0, 



NSPP/TSEP = 3.5?3. 


WTP s 738,0, 110.0$ 

•• 


SNAMV 

LSM = 2.38, INPRT = 1500.0, CNA = 15.47$ 



IPT TH 

STAGE 1 THRUSTING 



5MAM2 

ARPA = 1.5, lOOT s 1, 




TIM = 0.0, 0.01, 0.04, 0.05, 

0.09, 



.15, 0.84, 1,14, 1.74, 

2.04, 



2.34, 2,49, 2.64, 2.79, 

2.99, 



3.09, 3.21, 3.28* 3.3^* 

3.40* 



3.46, $,54, 




NTh/THS = 1023., 26439., 39417., 

42$97., 



42580., 43263., 43551., 44705. » 

45667., 



46820., 46723., '45570., 42974,, 

37975., 



33648., 24034., 1634$., 1086*4., 

6538., 


326P.» 0.$ 

fMAM3 

TWT = 0.0, 0.01, 0.04, 0.05, 

or. 09, 



.15, 0.84, 1.14, 1.74, 

2.04, 



2,34, 2.49, 2.64, 2.70, 

2.99, 



3.09, $.21, 3.28, 3.34, 

3.40, 


3.46* 


NWT/WPR = 738.0, 737.26, 732.10, 729.98, 

721.20* 



70$, 87, 553,07, 455,0?, 346.63, 

275.78* ^ 



203.27, 166.60, 130.42, 95.72, 

53.41* 



34.69, 16.60, 9.22, 4.95* 

2.22* 



•ftSf #00$ 



«mAM4 

SAREA r 1.4741, 




MCH = 0*. _ 0,75* 1.00* 1,20* 1.60* 

2.0T’ 2.40*' '2.80* 3.30** 4.00* 



■ 


NCD/Cno * .675. 


.925, 

.870, 

_ .780, 



.710. 

.660, 

.615, 

.565, 

.520, 


A 

.A55« 






CnAST 

.stage 2. CO a ST TNG 


• 


■" ' ■ 

<RmAM4 

-5AOPA s_.2?» 







VCH * 1.0. 

_2..p^ 

3.0, ■ 

0 

5.0, 



6.0. 

8.0, 






_NCD/CDD s .93, 

.72. 

.59, 

.51, 

.45, 



• A2, 

.36$ 




■ ■ 

2M0 TH 

STAGE 2 thrusting 





5WAV2 

FRFA = .223, lOPT = 1, 






TIM s 0.0, 

0.04, 

0.08, 

0.15, 

0.50, 



.70, 

1.00, 

1.20, 

1.40, 

1 • 60 1 



2.00, 

2.40, 

2.80, 

2.98, 

3,04, 



3.16, 

3.26, 

3.40, 

3,50V ^ 

5.00, 



nth/ths s 6., 

2050., ■ 

7400., 

7200., 

7075., 



7200., 

T700., ■ 

8080,, 

8275., 

8400., 



8500,, 

«7?5., 

9370., 

9050,,' 

6900., 



2700., 

oOO. ♦ 

225., 

50., 

0.0$ 

1 1 - 1.- 

$WAM3 

_IWT_. =. 0.0, 

6.04, 

0.08, 

■ 0,15, 

0.50, . 



.70, 

1 • 00 f 
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^royMOA*>STdP^/TWT*NWT»WPP ♦W^fdHTVMAS'^ 7 

CnMMOM/STOP6/oO«WPM,woL*WOo? *WTM 

CO MMOm/STOPT/y )0T,YnnT,Z'^0T*XnnAT,Yr>nAT»ZPOOT 

CdM V0V/5TdPA/riMPr:A ,PF,0A , pp ♦no TriM , a'Qm “ — 

CPMmOW/STOOQ>S,MX “ 

f O mmom75 too t 0 /s t 

T/'fHCn'MTT^^^UTrTTHxTTH 

^mmOn/STOPI 2>Yl.»VvL,'DYY*nAMLjPLnn ' ' 
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C0MM 0N/ST0P16/SA»APRAY>NN ■ 

c ommon /ST OR! 7/LFNnTH»<PL ♦OTSTtLVPL ♦LACC 

common / StOP18/TFPAr»LTNP* PT^T mp , CON V « I pNO » T ®P t NT i NP AOP ^ if P AG 

__common/storip/ttmp,to ■ " - - - - 
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C0MM6"N/SfdR"2‘4/PHTVNln,NfAB,NPH,KPH 
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COM MON/STOP26/ KCON»KSPENT _ _ “ 
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J_ NDLCHfZDLCH _ _ ' 

Cb>M0N7'ST0R2P7RR“An7AZPAh;FLRAD,VPLPAd;pRAr7NRAo;zRAdV‘^nPAnV' ■ 

J. NOPAOfZOPAO 

COMMON/STOPAO/DFLTA»VFL,ACrT ' „ _ _ ' 

__CdMM0N/ST0R81/SPHn ,COHll ♦LAMOAi ' 

_COm'mON/STOP32/NSYS,COV 

_C6mMON/STOR?3/NT IP,NCYC*NLTN “ " 

_C0MM0N/ST0R34/KLrN ~ 

'C0MM0N/ST0P3S/NTAP 

_COMMON/STOR36/PRESO*nFNSO*TPMPO*SOUNOO.Kl »C0M1 >KPRFV 

CdMMON/STOR37/TSLOPF»HALTB»TFMPB»PPFSR 

r'FORMAT( 20X.5H EL **F7,3,SX,5H AZ =»F7.3.SX,10H PAYLOAD s, 

i F7,l//I . 

I FORMAT (12'A6) “ 

! FORMAT (1 Ml) " - . . 


I 




i3-FOP'WAT(5X*12AA//) _ ' _ _ ' 

^dtmpmstom _ 

DTMPNSIOM TSLOPFC 0)*HALTB( 9J,TFMPR( 9),PRFRB( 9) ' 

OiMPNsrOM SMAT(3.^'i fTMAT(3*‘^i *S(o.6J ,5;T(2n.inj ^ 

DtMPNSIOM SoS(15»5) “ ' 

^DImpNSIOM_SA(^ 0,^10) ,ARPAY(10*30*lor‘ ' 

Ot MFNSIOM PHT { IT) *NIr>'{ 1A ) *WTM( 6)TTSEPT6 )T».n'TPTM 

^OIMFNSK^^TIMYSO) .THSj^SO) »TWT(50 ) *WPR( 50 ) *MCH( 50 ) *CnD( 5oT 

^bYMPNSTOM mfssTbO)' “ 

^MFNSl ON NPA'OE ( 1 O') — 

01 MANSION CON { 10*7) 

TNTP6FR THCON' 

AL J L Ajr, r t PLON , T T^, HPT I^T~ ^ 

A^^ L AMr)A»MAS S» LONO * L AT » MACH t MCH f L A vjS A 0 ,XAM 0 A I VLB AR ' 

RFAl K1 - - ^ 

real length ■ ~~ — 

PfAL_LSM»_I^FRT' ' 

RFAL NR'AO , NLChVNDR a'O • NOL CH ~ 

DOUBLE PRECIS TON CC'S" 

EXTFRNA'L ACCEL1*ACCEL0,ACCFL3’' 

_>^AMO/VT,FL,AZ,NGFO*PHr,LAMOAVTIM{;,A(.T,LFNr,TH ' ' 

name*. 1ST /NAMR /P T T ; NSOPNT , NS TOP , T STOP , NSK I P i X STOP~» 

/ 1 tH rONttl»'ALTC,NSY5;»NT IP^NPAGF " 

! NAMFLT ST /MAMS7PHT0»LAMnA0, _ PHT_1 .LAwnAf 

_J NAM<^LTST /NAMy/NPH5*PMT,NT4R*NTn 

N AMFL l ST '7 n AMU/NMT , WTM , VI PL , NSEP « T SPP t WT P 

NAM^ Yt5T/NAMV/L5M, INFOT^CNA - 

NAMFLYST /NA'M5/>Ti VWTF*NTABtNiO 

pAjA BLKK/b20?0202020?0/ “ — ' 

data A P 0/6HAP0(^ EE/»TST/6HTST0)^/» ALTO /6HALT=0 ~ / " 


c tnttialtzf constants 

p L..= . 3 . 1 4 1 5 9 2fi 5 

OMFrA = 7,292nF-A 

CON V = PT/ 1 e 0 , 

P T 2~ g 2 « »P T 

TO = .0 " _ri; T 

P0'' =~2U5.666 

EPTfNY' = 5.0F-7j;3_ _ 

EPBTG = T.0E-i5_ 

LOW_f 1*0E-4 

GM = 1 • 40765 7AFf6 

RE g 20899262. 
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_RA ■»'''209?57A1.' 

■ ■* ” RB a POS'^SSOI, 

'\6p s ‘?2.174 

* ALTM = 2.0P+'' 

LTNP a 0 

IFNO a 0 

Ncvr a'o 

NLIN a 0 

KPAA a 0" 
ifCON a T" 

KSPFNT a 0 

NAP" a - 

rspT a 0 
bo 14T Jal,14 

% ^PHT ( J) a' 900000. 

141 CONTTNUF 

DO 142 JaUlO ■ 
NP'a 7;E(J) a 8~ 

142 CONTINUE 
■pHAS'a BLK< 

DO 14' I=l*lO 

bo i4 j=r*7“ 

( CON(I*j) a 1.0 

140 CONTINUE 

, CONT5*l)'"i .3048' " 

jCONiS.l) a" 1.852 

CONI 1 >2) _a_4.4482 

CON (2 ♦2)'=' 4.4482 

CONI3.2) a 4.4482 
CON(5.2ra 4*».88 

CONr6»'2T'a ".3048" 

CON(7,2) a .3048 

C ON(B«2) a 14.594 
~~C bNri»3r=' .3048 

CON (2, 3) a .3048 

^CON(3,3)_= .3048 _ 

CON (4 »3)' a “,'3048 

CONI 5 *3) a .3048 
C0NI6*3) a .3048 
CONI 7*'3 )' a ”.3048 ' 
C0NI9,3) a ,3048 
CONll »3) a .3048 
CbNIl»4) S”r.852 ~ 
C0NI4*4) a .3048 


I 




I 


( 


c 


c 


c 


_CON(5,4) = 
C0N(6*4) = 
-CON (7,4) = 
_CON(8V4) 
CON (9, 4) 


1.P52 
1.852 
J.852 
= .^048' 
= .3048 


_C0NJ[1_,4) V ,^048 
CON f 3 , 5 ) '^’17852 ■ 
_COnV1 *6 ) _= _. 3048 J 
CON ( 3 ,6 ) g' .3048 ~ 
CON (‘4, 6 T = .3048“ 


REAP OEN'PR AL ~ T NpQT DA t'A FROM "T ARP 


INPUT 
-PRINT”! 2 _ _ ' " 

REAP (5, 11 r fi Tie 


R PAP (5.NAMQ) 

RRAn(5,NaMpr 


REAb( 5,'NAMS) 
_RFA_D(5',NW)I 
JREAD(5,NAMU) " 

_RiAn'(^'NAW> 
-PBINT_13’» title “ 
print 10,El,A>,WPL 

Sf PTI = FT I 


>00 


_JF(NITP .FQ. 2) NTftP=>~ 
__IP(NnP .EO. 0 ,0R, NT TP 
NCYr g 1 0.»PTT 

-p-iT = .1 r“i 

^CONTINUE 

JLIO^* PTI/2. ' 

_pT0 = PTT>'8.— 

0t0_= AMTNI (0,1,DT0) 

DT g ' f)T0 __ _ 

-Ik -®- el*con~v ■ ' ■ 

AZ*C0NV 

-PH lQ ° PHT O»C ONV 
EPS = PHI 0 " 

-P^ij * PHn*CONv' 

PHI = PHT*CONV 


.FO 


. 2) '^O TO 200' 


LAMOA _«_LAMnfl0#C0NV~ 
L AMOAl _s LAMOftl^cOwy 
LAmfA a LAMnft*rONV 


CA2 = 
SA2'» 


COSIAZT 

SIN(AZ) 


(. 


c 

C c r^EODPTIC TO OFOrPMTPIC rONVFRSTONS 

APO s( RB/Rft ) **?*s I KM PM? 0 ) /ros I PH T 0 ) 


f fPS.=_fPS - PHTO 

AP Q g t p R / P A ) **2 «lS T N ( OH I_ ) /COS ( PH U 

PH T = A T A N ( A RO ) , ' 

C AP O_a (RQ/PA)»»2»STN(OMn J/rOSfPHTlV 

PHTl g ATAN(APO( - 

SEL » SIN(FL) ~ 

C SELGD « SEL " ' 

CEL = C~OS (FL( ' 


C C EPS = C OS ( EPS ) ' ^ " 

SEPS s S TNfFPS) ' " 

SFL = SEL»CPPS - SepS*CPL*rAZ 

C PiL^ATAM2(SPL.SQPT(l,0 - SPL**2J ) ‘ _ 

AZ_ g ATAN2(C FL*SAZ «C5 pS*CP| »rA> + SFPS*SFLG0) 

# A ft M ** ** ' I I mm^ I I • ■■ I Ml ^ 


SAZ =' SIM(AZ) ' ■ 

?so con tinue ' ' 

c / c:zi__. “ 

’ c ROCKFT MOTOR WEIGHT 

TF(NMT »fq« iT'go to ^*16 

C SUM f _ .0 " " 


SUM = SUM + WTM(J) 

C rLP-Co'njJnuf 

WPM = SUM 

__NMOJ = 'l ' 


C TNltlALTZf WfnOMT f'P PP00FLL*MT 

SUM a ,0 - - 

( DO ... .31 5 J= 1 *NMT’ 


RlS^CONTlNUE “ 

( ”_WPP.?. r SUM 


1.6. WPP 2 =' WTPrf) 

( ^’WRM'j* WTM(l) 

~N MOT i"' _'l 

-a 17 CONTINUE 

( ' IkDFL a 1 


I 


( 


NOPL = 1 _ 

( KPH = 1 

NPH = 1 _ _ 

c " ■ 

( ■: READ FIRSt PHASP D&TA 

CALL TABLF 
Npif s NPH + NnFL' 

C KPH *"KPH + KnEL 

SPHfO =SIN(PHI0) 

CPHiO = COS ( PHI 0) 

( ~ SPHtl* SIN( oHIl) 

CPHTl * COS (PH ID 

XO « RP*cWfO*COSaAMnAO) 

( YO s' RE*CPH lO*SIN(LAMnAO) 

20 = RE*SPHI0 

SL s SINTLAMDA) 

( CL s COS(LAMDA) ' 

SPHT = SIN(PHI) 

CPHI = C0S(PHD 

( LONG = LAMOAO 

^ call launch 

\ Vx = vt*co$(El‘)*co5;(A7)’" 

( _ _ VY = VT*C0S(FL»*SIN(A7) 

i VZ g~VT »SIN ( el ) 2 U 

i : IFINSKip .E0,“ Oi 'Gn' T o ' ife ' 

( R = RF -f ALT 

LB^R = LAHDA +^MEGA*(TIMF j^TO) 

Z = ’R*SPHf 

( 1 J = R»CPHI*SIN(LBA0) 

2 X_a J?»CPHI*C0S(LBAO) ' 

_CALL HATROC " 

( 1 NLTN = -1 

GO TO _ ■ 

^^18_C0NTInUF‘ 

( 1 _CAZ = COStAZ) 

§AZ_*_SIN(AZ)_ 

“VXL = VX»CAZ +VY»SAZ 

V C ■„ _ VZL * VZ 

f X L « length»cfl ' 

i 


( ( 


ZL = ALT 

theta s cl 
call' ATMSPH 


2 7 SOLVE AN6LF OF ATTACK EOUATIOMS 




0 


’if 


M 



( 




i: .. 


c 


( 


2 a U 2 

( R_e^jhRT(X**2 +_Y**2 + Z»*2)_ 

ALT = P - PF 

VX = V XL»CA2 

( V Y a V XL»SAZ 1 

VZ a " VZ L ■ 

<<10 CONTIMUF 

c -C-Z z. Z 

C COMPUTE INTIIAL VELOCITY 

CALL t R AMS‘2 ( VX . V V* VZ * 1 »X ♦ U Y ♦ U'Z ) 

( XOOt'a^-nMEfiA^Y + UX 

YDO t = 0MEGA»X » UY __ZZ 

^?noT uz ” . _ 

C TNfflAL TZE s" ARP AY ' ' 

N a 1 

_STi/N) a X 

S(2>N) a Y 

SC3.N) a Z 

S(4*N) a XDOT 

S(5 »N) a yDOT _ 

1 S(6*N) a'zoof 

c (C ompute accelpratton _ _ 

_NX a“i 

CALL ACCEL3 

LHETA a PL „ 

ELPOFV a FL 

2 ' ^^VT ,LT. l,0E-5)"6f> TO 320 " 

c‘“ ' ■ “■ ' * - 

‘C COMPUTE OUTPUT DATA 

■ CALL OATAl 
' NN a 2 

c 

C STORE DATA TEMPORARILY 




( 

( 

C 

c 

( 

i i 

I 

t 





CALL A^LFA . . . 

CALL MATLCH 

Nl^CH ^ X L*CA Z 

ELCH a XL*SA2 

2LCH = ZL _ . _ 

Z ~Z a Z LCH -_RE ... 

CAL L _T R AMS 2 ("nl'cH * EL c'h , 2 2 * UX ♦ U Y , UZ ) 

X a UX 

Y a Uy 


■ I 




c 


( 


-CALL-SETSA 

_IF-tNS<,lP_.EO. _0) T’'-> ^20 

PTTM E 8 T T_Mf 

Xffac =.j ?*0 

J^.UP.STl=rN.ITp.’. ' Z __ 1 

NTT P = 0 

CALL out 


NI TP g N TIPST ^ 

TTMi a ATNTCTTMF/STPTt 1»SfpTT + STP'fl 


TT M? a A TN T(TTMF /otT ) *PT I PT I 


- t~!M 2 * l>Qg~ 5 n"^~ 


*^20 CONTINUE 


PTTME = A T»JT(TTM P/PTTT*PT t + PT T 


TPREV = TTMP 


N = 1 
ST(1*N) 


C 


I- 


c 


a X 


STC2*N) = Y 


STO.N) = Z 


ST(4*N) a XOOT 


c 


ST(3*N) a YDOT 

/ 


ST( 6 *N) a ZOOT 

\ 


ST( 7 ,M) s XODot 

c 


ST( 8 tN) a YDDOT 

7 

ST(P*N) a znor>T 



^IF(VT ,GT. 1.6'P. 

( 


OT = LOW 


call rungeuccfi 



OT a otO 

( 


GO TO 351 

c 



c 

FI NO “present' 


r.O TO ^^0 






CONTINUE 
JFjtABSd.O - 
GO TO 332~1. 

Wrm“= WRM WTM(Nynit j 

-SPT a LSPT 1 

1**> 

SPs VJ«LSpT~) V ST uVi ) 


. ROCKPT MOTOR WEIGHT 
TSEP(MM_0 T_)/TIME)_.L^^^^^^ go to 331 


- ( 


^30 CONTINUE 


_SPS(ll.LSPT)_r el 
JPSd tLSPT) s TIMP 
NMOT e 'NMOt" + 1 


332 CONTINUE 




c ' REAf) TN data” for Ne'xf PHASF_ ... 

_ _caLl table 

KPH_f KPH_+ KPFL 

“■ NPH = NPM + ’npFL . . 

_ CO NTINUE _ _ 

c " ^ '* _ _ 

c talculatf next TRAJ'^CT0RY_P0.TNT at T 4 dt 

' call' PUN rnF (ACrFLl ) _ _ 

•^51 C ONTINUE _ 

TPOT '= t'f ME DT __ _ _ _ 

IFitTPOT-PTIMF) .6T. LOW .AND* tiMF -LT. PTIM^) d T s 

1 PTTMg TTMg'V LOW -- - - 

C c compute ~ow . 

CALL OAT A1 

; NN a 1 ~ ' ■ 

( ^ STORR ■OAtA'trMPORAPILY ~ 7 ”1 

' call sftsa' ■ 

( • - "ipj^STOP .EO. 0)' r,0 TO 708 “ * _ 

; ■ r 'v. • 

* C CHECK STOP TIME 

I <; IFITIMF .OT. TSTOP} r,o TO 703 * ” _ T Z I 

f- iF CABSn.O ~ TIME/TSTOP) ,LT. 1.0F-6I G O tb~70l ; 

• _g 6 to 70fl ' 

■ ( _ Jl’jOl TENO a 1 : 

T FP AC a I ,0 i 

UPTIME a TIMP ’ ~ . - . 

60 TO 801 ^ 

‘ tOVlENO a I ‘ _ __ _ ; 

TFRAC“a '(TSTOP TPRFV)/(TTMP TDRFVI ' 

( _ NTIP = 0 

- PTiMP _ TSTOP _■ J" ■ I 

> _ PHAS a TST ■ ■ ■ ■ ? 

( _I _ 60 TO eOl - . - - - ^ 

•rOe CONTINUE I 

~C _ _ ’ ' ' I 

i ( ^ ' ALTlTUDFaO CHECK “ _ I 

i .AMO, FL .LT, 0.0) 60 TO 710 ‘ | 

_ — _ I 

f^( ■ 7~lO'TFRAC a SA(S,?*l)/(SA(S,?,l) - SA(A,1,I n | 

'"lENO a 1 . ‘ 

f >>T'IME a'TPRFV + TFRACMTIMP - TPREV) 

( NUP a 0 




I 


i PHAS » alto . 

i . i30 ...TO 80l 

7 ?0 C ONTINUF . 

If} MAP .PO. 2) r-0 TO 

( 

C A Ppo F P_CH 

^rF(FLPRFV .Cl,_0j^0 . 00 . -fL »G~T. 0.0) r,Q TO 7?1 

( ’ PHAS = ADO ” ~ ' ' 


?c 



KOFL = 

1 


NDFL = 

0 

IFTKSTOP .EO. 0) GO TO 728 


lENO = 

1 


TFRAC 

= 1*0 


PTIME 

= TIME 


^ GO TO 

80l 

728 continue 


NAD = 

1 


tfpac 

= T.O 


PTTME 

a TIMP 


NTTPST 

= NITP 


nitp = 

0 


CALL OUT _ _ __ 

n.itp„=...nti.pst__ . 

A.T NT IT.T MP/S TPT T. ) *S.TPT. 

_TIM2 = ATNT(TTMF/PTI )*PT I +. 

NLTN = NCYC,.-. TNT(lO,#(TIMi 

PT Tf P _ = AT NTCT tM F / p T T) *PJ T_Ji 

GOJTO 80s 

._t2.1_C0NTINUF 


I _+ 
PTT 
- T 

pt t 


.STPJJ 

I M2 + *lTOf-SjT. 


TF NF'xT PHASF has FTAPTFO, PRINTOUT OATA 
TFflNFXT ,NF, n 00 Tn 7fiO 

INFXT r _0 _ _ ■ ■ 

PTT ME » TfMD “ 

_ __CALL OUT 

' NTTP = NIIPST __ ~~~ 

^flMl r ATNT(TTMFVstPTt )*STPfl -i- STPTI 

JIM2 = ATNT(T1MF/PTI )*PTI + PTl 

2 _NljN_= NC YC - I NT(lO. »(TTM j ~ T IM2 -jh"! .o'p- « 

DT = OTO - . - - ^ - 

"PTTME » AINT(TIMF/pTn*PTI _+ PTT 

^g6_to 80s ' 

780 CONTINUF' 


( 



( 


AT PNO OF PHASF* POTNTCIJT DATA . . . .. 

_TF(TNPXT ,NP. 2 ) TO 7P0 

IF( A BS( 1.0 - phT (.NRH) /TJMQ_,.(\T, l.QF.S ) 0 0 TO 7Q0 

TFRAC .« 1.0 . 

PT.IMF ss TIMF _ 

NTTPST » N I.I.P 

NLU».».0. 

CALL_0UT. ' 

LNFXT_*_ft. 

PT LME „• A T NT ( T 1 MF /PlU *PT T._+_ PJ.T 

00 JO BO'S 

79 0 COMTINUF 

~ rHP^ tVmf Vb~R 'pR ui~fnuf ' 

! F(TE N D_. EO_. _l.L 0.0 _TjD_ .8 0 1 

^IFdlME .LT. PTIMF) C-o TO 805 

IFfTl ME .FO. PTIMP) r-O_T0 800 


COMPUTE TNTFRPOLATtON FRACTION 

^TFRAC_a (PTTMP - TRRPW)/(TTMF -_TPPPV) 

GO TO 80 1 

.? 90 _TEPAC..a. 1 .0 

PtiMF a TIMF 

rOi continue ' " LI’ 


CALL PRINTOUT ROUTINE 

mu_ouT 

IFUENO .EO. 1) GO TO 840 


SET NFW PRINT TtMp 

PTIME 'a ptYmE +'PTI 

R0*i CONTINUE 


STORE PRCSPNT VALUES OF ALL OUTPUT DATA. 
_ NN a 2 

call .SFTSA . 

__TPRFV a TIMF 
ELPREV a FL 

8 06 CONTINUE 


PHASF TIMF CHFCX 
' ~ IF ( NAP .PQ. 1) GO TO_R10_ 

go’ TO 81 V 

ftlO NAP a 2 



i 


DT e 10,0»LOW ... 

( GO TO (340,332)*KCr»N ... 

Bll CONTINUE 

~ ' IF(NPH .GT, NTAB) r,0 TO 330 

( 1'1TF(ABS(1.0 - PHT(NPH>/TIMF) ,GT, l.OE-5 )' GO TO 723 

KDFL e 1 _ _ 

_ NOFL =1 _ • 

( INFXT » 1 . _ 

DT jr 2.*L0W 

GO'fO 340 ' ■" ■ " 

C -^2 3 'CONTINUE'" ' 

IF((T!ME >»• DT) .GE. PMT ( NPH) ) 'GO TO 730 

IF(ABS('i.O_- PHTrNPH)/rTIMF>~DTT) .LT.X'O^-^' GO fojj'fo 

( _ 2 _G0'T0' 350 “2' ■ * - 

*730 CONTINUE * 

IFTpHTINPH) ,LE. PfiMe)"‘GO“‘fO fiT 

( IF(TPhT(NPH)' - PTTMF) .LT« 4>»l6W) go ~T0~713~ 

D'f = PTIMF'-'tIMF + L<^W ’ ' 

GO To 3'50 _ ~ ~ _ _ _ 

C ^ 1.1.1 continue ■ - - 

0T '= PHT( NPH) '- time 

( INFXt ='2 ' ~ ‘ 

C GO TO 350 



; C ■ C 'spent STAGF ROUTINE 

B40 'N PH S g 0 

NP*H » i " ~ - ■ 

(, ^ kCON = 2 

' "TF{NSPE_NT)PQ9*999*R50 

350 KSPPNT e 'KSPFNT + ' 

( KPH » 1 

<OFL_= 1 _ 

_ KPAr, e 0 

( ~ DO 65 J » 1.T4 

PHT(J) a 900000. 

630~CONTINUE' 

( ~ _ 1" 00 «55 J a lf9 

i ' I F(NPAGE(J) .FQ. 5) GO TO B56 

.! R33 C'CNTINUE 

V ' ( _ " GO' TO 85R 

T'*' Jf»5'6 DO «57 l a j,q 

_ "* NPAGEn'*)' a ’N0AGF( T+Y) 

i { R57 CONTINUE 




I 


i 

I 

f 

t 


I 



C&rCPLO SUP. ACCFI.6. — - ALFA A^CFL.^ AT lON.FOS^, 

SUROOUT I wE_ACCFI.O 

boURLF PPFCisiOM' T?-*AT^SMAT 

_ C OMMON / S TOR 1 / A L T ♦ T F M P , PR E S 7d E N S_* V I S C » Sj?UN 0 f * 

C0mm0M/ST0R2/X»Y.7.»R._ 1 ] ' ' IL_J 1 

_C0MM0N/ST0R3/TMAT».SMAT»SPH1 *CPHI »CL».SL 

COMMON/STORa/SAPEA .MCH.CODjCn.MACHtNCn.CDS.iSRAG ^ 

C OMVON / S TOR 5 / T V T ♦ vop . WE I_GH T . MASS 

./ycoMMON/sfbR6/p6rwpMrweL»VpP2 .wirM 

COMMON/STOR7/XDOT .YDOT ,zbOT ,XDOOf ♦YDDOL»ZDDOT 

_ C0MM0N/ST0R8/0^'F^5A,Pc,RA»RB♦G0*'G^i»ALI^ : 

COmmOi^/STORP/S.^JK ^ 1 ^ 

^ ‘JCOMvON/STOR13/"xlVzL .THFTA .VXL_.VZL_.ALPH .'GAHA.CA2.SAZ ~ 

common / S T or 1 a / L ,J_M cpt VcnA 

C OMMON /sf OR 15/0. VX , VY , VZ_. VJJ EL_. ^ . RANG E . l_A_T .LO NG. ACC tGDLAT 

„_Z C0mm0N/ST0R1P/TIMF,T0 . 

’ _ COMMON/STOR20/ARFA. lODT.tlM.NTH.'THS.THPUSt 

1_''~ COMMOM/STOR23/SPHIO.COHl6.LAMDAn.XO’.Yb.Zfr” 

COMMON/STORZS/LOW.EPTTNY.EPBIG.PTT.PI .PI2.bt .HIGH ^ 

DIMENSION TIH(50) . THS '( 50) .TWT ( .WPR~( 50) .MCH ( 5Q j .CPD ("5 q') 

DIMENSION SMAr(3.3) .T'mATTB.S r.^M(j6jiSj^ 

RFAL LAT .LONG.LAMnAO.LOV' ' 

I_1_1’’‘^AL LSM.INFPT.MASS.M aCH~.MCH 

^ _ N_=__NK_2_^ 1 

’ ' XL * S(l.N) 1__ 

zi,_*_S.!ljLNi__ ' 

ThFTA = S(3.N)_ - 

VXL = S(A.N)_ __ 

yZL _S( 5 .N) f " 

. VL .8, S0RTIVXL**2_ + . VZL**ij 

CALL atmsph.. : 

MACH_ = , V.L/sbuHO 

...0 = (bFNS*VL»*2 )./2, 

C COMPUTE DRAG 

CALL TABI 

_ DRAG s cns*6 

"_c ' ‘ compute'thrust ^ ^ ' ^ ' 

' C Al \ TAB ? ' 

C . COMPUTE WFIOHJ 

CALL TAB3 _ ^ 

MASS a WFIGHT/GO 




1 

} 

i 

I 

! 



h 


\ 


GAMA = ATAN2(VZL»VXL) 

■ 

ALPH = THETA - GAMA ' 

CS = COS(THETA) 

. 

SN = SIN(THETA) 


FNORM = CNA«SAPPAw6*ALPH 

- 

AXL .=. (THRt'ST*CS - nRAG*CS - T- n6pM*«;M ) /MASS 


AZL = (THRUST<-SM - D9AG*SN + FNOPM*CS ) /MASS - GO 


thdo = (-l'sm*fn6pm)/imfrt 


S'(7,N) = AXL 


S(8,N) = AZL 

S(9,N) = THDO 


RETURN 

END 




) 

j 

'<f 




m 


I 


N 


riMO J 



C^^L'i _ SUB. ACCF Ll ACC>L>RATf 6 n' FOR RUNGF 

~ SUBPOUT IMP ACCFLl >: ; 

"noURLF PPPCTSTHM TvAT,SMAT 
c OMMOn / ST or i / a L t ♦ T FMP ♦ PR E S • f)EN S . V I SC ♦ SOUND 
COMMON/STOR2/'X,Y.r!rR' “ 

_cbMMON/STOR3/TMAT»cMAT;*SPHT .CPHT *CL*SL ~ 

” COMVOM /STOR 4/ S’a P P aVV!Ch *000 ♦ CD ♦ MACH « NCdVcDS«DRAG ' 

C0MM0K!/ST0R5/TV.'T,N'.-.'T»«*.pR,wEI6HT»MASS 

COMMOM/STOR6/PO»V.'RV,WPL»WPP2 *WTM 
C0MM0M/ST0R7>DXDt^^ , XDDOT . YDDO'T « ZDDOT 

■■■ COMVOn/STOR8/OME6A,RF,RA.RB,GO»6M,ALIM ” 

" C0M'^0N/ST0P9/S«NK ' ' 

■ COMMON/STORll/tHCON»tl VALTC»THX*THY.THZ 

COMMON /STOP 1 5 /Q . Vx7 VY » VZ » V^* EL . AJ, RANGE .1 .LONG , ACCtG D LAT 

d; ' COMmom/STORIP/TIMc.TO 

common /STOR 2’07aPE A , lOP fTt IM ♦ NTH.T HS .THRUST 

~ COMMON/ STOP 26/ kcON ♦ K FPFNT* 

^ DI mfNS I ON TIM (' 50 ) ,YhS fsOJ VTWT ( 50 ) ♦ WP^( 50 ) ♦ MCH ( 50 )'» ODD ( 50j »W fM(6 

^ D 1 MFNS I ON SM atY 3 ^ i »TmaT( 3,3) »~S(9«6) 

'INTPGFR THCON - - 

' 'real' LAMOAO, MASS ♦ LOW, LONG* LAT, MACH, MCH 
N = NX - 

INPUT present' POSITION AND VELOCITY 

X = S ( 1 , N ) ~ ■ 

Y = S{2^N) ' 

Z "r Sj- 3 t I - 

bXDT = S(4,N) ] [ : 

pYOf ='S'('5,N) - 1 

CXDT = dxdt + 6mfga*y ^ ^ 

CYDT = DYDT - OMFGA^X__] ^ 

DZOT , = S ( 6 ,N ) ^ . 

CALL MAT ROC 

CM. L TRA NSI (CXDT, CYDT, DZDT,WX,W Y,WZ) 

yx_.?_WX CvgLo&t^2A^ 

VY_=.V/.Y ; 

yz_ = .wz : 

yi_=_ soRT(yxf*?..+ vY *»2 .+_.yz** 2 .) ’ 

R . ..SoRT(_X**2^+ Y»#2 i 

ALT.= R_- RF ___Z 1_ ^ • ‘ 

IFCALT .F-F, ALIM)_ GO TO 80 _ _ _1 

CALL_ATMSPH_T02FJNn S^ED'OF"' SOUND,' DENSITY,' AND^PRESSU^E I 




CALL ATM^PH ^ 

MA_CH = VT/SQUND ;; 

...1 GO ,fO0 ! ^ 

80 . CONT i NUE ^ L__ 

MACH r_-2.0 , 

PPFS = 0,6 ; 

OF NS ~ ='..P , p ' 

■ • ■ 2 

Q_ = _ In P N S* V T .♦ ♦ 2J_/ 2. • •• 

C roMpy.TF_r'PArr 

C ~ rOMPUTF .thpusj 

C COMPUTE TOTAI. WFIOhT 

GO. t 6_ ( 9«> /gs )_♦ KCON 

- 55_CALL _TABJ^ ; 

** CAlL TAB 2 ^ 

CALL TAB-^ ■ 

GO TO* lOO 

9.8 CALL. T A B ] ^ ^ 

Sf . 0 , 0 ^ ; 

WFTr,HT g i/TmT<BPPN^ ^ 2 

.1^6 Onwf j.MU.F ^ 

' 1.’ DRAG = COS*o’ ^ .1 

f MASS = WFIGHT/GO _ 

' C_ COMPUTE GRAViTATiONiAL ACCELERATION ■ ^ 

6G..= _GM/P_»«3 ^ ~ 

.rGG*X 

6 Y ..= - 6G * Y ^ 

6 i =. rGG* z __7 : ~ __L_1 L 

FX. = DRAG*yx/VT 1 

.FY = nRAG*VY/vf I 2 l 

FZ'= drag#vz/vt ' ’ 

CA LL J.P A M S 2 ( FX « > yVfZ «"xThY ♦UZ ) ' <l'r< KO " ~ 

i CALL DIRPCT. _ . 

C COMPUTE TOTAL INEPTIAL ACCELERATION C0Mp6n'ENT.S 

XDDOT = GX 4 ( THX - UY ) /MASS ' ' 


Y 0 DO T _=_' g Y_ 4 '_ ( T H Y ' Ju Y / >yM A S S ^ 


S(7,N)_^XOpOI _! 

S('8»N) = YDDOT _ • • ‘ 

i S(9,N) = ZDOOT ~ ' " • -- - - - - - -| 

' " . RETURN : 7 7_ 77. J 7_771 

■l.-^ n .d ~ ~ ' r i 

?-v • ' 



X 

4 

o 


.CACCEL^ JS-UB._ArCF.L2_r^==^LAmCH_RAJL, ACCPLFRAJIQN „ 

SUBPOUT I NP ACCP-L2 ^ 

ppURLC_PPFCTSTON ST 

COyMbN/STOPT/AL'f ♦ ♦.PPES ,DFNS ISC ♦ SOUND ■ 

C OMMON /ST CP 4 /s & R C' A , f/Qu , CD D » C D » M Ar H VnCH ♦ CDS » DJ? AG 

COMMOm/STOR5/TWT»NWT*J».'PP*WEIGHT*MASS 

'COMMOM/STORS/OMEGA »RFrR'A »rb,go\gmValini ^ 

COMMOM/S’TO'’Q/SVnK 

M 0 M / s T OP 1 o7s T 

COmCio'n/STGRT5/0,VX.VY,V2,VT,EL,AZ»RANGE.LA T.I^NGj 

COMYON/ST6Rl7/LENGTH*SJL_»6lST*LVrlVLACC' 

COMWON/STOR19/TIME»TO 

COMMON / ST OR 2 0 / A o_E A V I OP T , T T M ♦ N TH * T HS , THRUST 

D^MPNSIOM ST(iO»lb) ^ __ 

MCH(bd) .c'nb(SO) ♦TWT(50'i ♦WPR( 50) ♦TfM('5oi 

^EAL LVEL*LACC*MACH*MaSS 

N_=_Nlf 1 

DIST = J 

Ly F_L , 5. _S. ( 4 » N )__ ~~ 

MACH.= LVFL/SOUN D 

Q = jbFNS«LVFL»»2 ) /2, 

C _ COMPUTE DRAG 

CALL.jABi..l’_; - 

DRAG__= CDS»0 

c computfIthrust " 

.CALL T A P.2_ ^ 

C rOM PU T F WEIGHT 

CALL.TA.Ri 

MASS - .W£.IGHI/GO 

.LACr - (THRUST_.-_WEXPMT*$FL_- DRAGr/MtoS 

...S(7,N) LACC ^ __I 

RETURN 

END ■ ^ 


9 


I 


I 


o.c. 


_CA_rCEL3_ SUB. 

1._SUBP0IITIMF 

pnUBLF _ PPFCTSIOM tmat.smat 


_ACCEL3_-~_ 

AfCFL3 


ACCFLFRATTON USING EL»AZ 


COMMON /STOR 1 / AL J ♦ T Ff^P PPFS.r>FNS»VlSC» SOU ND 
COMMON/STOR2/X*Y*7~,P ^ ’ 'll 


C0mm0N/ST0R3/TMAT,5MaT»SPHI*CPHI.CL*SL_ 


COMMON /ST0R4/SARFa,Mch, CDD»CD» MACH .NC n .CDS ♦DPAG^ 
COMMON / S T OR 5 / T '.n » N V T,» W PR *WE_IGHT,MASS 


. C OMMON /.S TOR 7_/ X DO T , YDQ T ♦ZD QT^ XOPQT »YDD QT»ZDDOT ~ 

_CO_Mm6n / S TOR 8 / OM F G Aj P.F, R A , R B »G_0_,1 gM A L LM 

..,C0MM0M/ST0RR/S*NIC 1 

common /STOR 15/0, VX , VY ,'vZ » VT rELVAZ', R AN g'E ♦ l AT , LONG » ACC » GDL AL 

COMMON / S TOR ? 0 / A R F A , I OPT , T I M ,LTH , T HS *jr HR UST 

dimension SMa't(3, 3.) ,TMAT(3,3) »S(9, 6) 


-P-I, Mf !1S 1 0!LJ I 50) ,THS( 50 ) , TWT( 5 0) > U’ PR( 50) »MCH( 50 ) .Cnp( 50) 
JRPA.L. _ M A S.S »j-PNG » L A~t, MAr H,MCH 


N K 

1 N PUL. PR E S FN L PPS IJJP N. A_N P V FLOriTY 

_ S( 1,NL 

. _Y. SA Z* N .)_ 

Z = S(3 ,J11 


{ 


PXOT = 

DYPT_=, 

CXPT. = 

CYD,T_ = 

OZPT_= 


S{_4,N) 

S{l5,N) 

OX n.T . . + , OM F A ♦ Y_ 

OY OT - 0MFGA*X 
S(6,N) 


.CALL,«M.RO.C 

CALL .T.RAN.SI (CXDT»CYJ)J:»D.ZPJ.^WX,_W.'L,WZL. 


-_-VX_=_WX_ 

^VY_ = 

__ VZ = 


WY 

WZ 


.VT..= S0RT{VX*»2. + VY*»2 4- V7**2) 

XELX^Z__t_CQS led^cosT az) 


CELSAZ =. C0S(EL1*SIM(AZ) 

SEL *.-SlN(ELl 

IFt VT. ,LT . 1 .OF-5 ) _GQ .JjD__7.0_ 

CFLOAZ = VX/VT 

.CEL5A7. = VY/VL 

SFL_*rVZ /yi . 

.... continue 

R_ =l..SpRLf_X **2 Y **2 Z»»2) 


70 


- 1 


n 

i 


( 


i 



_._.ALT =.R - RP : 

_L£.(AL_T„, (IE.. /LLXM )_ XLO_L0_Q .0 1 

call ATMSPH to FINr> SPEE D OF SO UND. DENSl TY ♦_JkllO„^ESSyRE 

CALL ATMSPH : i 

MACH - = - V T / soy Nt) 

60. TO 90 ^ 

80_CONTlNyE • 

MAC H = »?.0 ^ • 

PR ES._*„0 *.0 ^ ; ^ 

.. -DENS. = 0,0 : 

90- CONTINUE 

.0. = .(nENS*VTJL*21_-^, 

-CALL- TAB! : ; 

jaP-Ao = cns«o ^ • 

CQMPUT e thrust 

XALL-IAB i ^ :: ’ . 

COMPUTE POC KET WEICHT 


call. -T-AB-?, 

.M AX5._= we I GH-TZIlO 

DTFF = THRU. S7 - OPAG 

COMPUT ^ gpavt tattom acc flf pat lqn 

6G = 6M/P*»3 : 

-6X. = .rG6*.X 

.6Y . = ~GG» Y 

_6Z_= . -G6»2 


COMPUTE A FRO FOPCFS TN TOPOCFNTRTC SYSTE M 

-JLX . .= . D I ff»cflcaz 

-FY = D1FF»CFLSAZ 

. -FZ. =-0IFF*SFL 

^T.RANSFORM AFRp_ f6p.CES._IQ INERTIAL SYSTEM 

...CALL TRAwS2(PX»PY,ez,nv,UY»U7) 

CO MP UTE INEP T I A| ACCELERATION 

_XDDof_ = _GX + UX/MAES 

YOnof 's GY V UY/MA.S S L__ 

_ ZDDOT =1 GZ“+‘ UZ/MA5S 

5(7 »N) = XDDJDT 


S(8,N) 
-.S(.OjNl. 


„ . re.turn. 

_.....EN0..„_ 


= YDDOT 
..^ZDfJXLL. 






CAl FA SUB. ALFA ANGLE OF ATTACK EQUATIONS 

SUBROUTINE ALFA 


POURLP PRECISION ST 




DOURLF PRECISION OPl 




COMMON/STORl/ALT*T'^MP.PRES*nENS*VISC*SOUND 



0 . 

C0MM0N/ST0R2/X,Y*Z.R 

- 



COMMON /S T0R4 / SAPE A , MCH ♦CDD t C-0.» MACH » NCD »CDS»DR AG 


..C OMMON / STOP 5./ T W T_, N W T , 5.'pp_, WE I.GHT . 

COMMON/STbRT/XDOT ,YpOT .ZOOT , XPOOT tYDD Q T >Z DDOT 

COMMON / STOR8 /omega , PFj RA ♦ RP.,.60.»0M*A L 

C OMMO.N / STOP 9 /Sf N K 

.COMMON/STORlO/^T 

.COMMON /STOR 13/ XL ♦ ZL » THE T A »VX L » VZL » ALPH «GAM A,CAZ tSAZ 

C OMMON / 5 T OR 1 4 /^S^.J N FP : 

COMMON/STO'R 1 5/0 1 VX , VY , VZ »\^ * EL » AZ > RAN GE »l AT « LONG> ACO.^pLAj[ 

C0MM0N/ST0R16/SA,ARRAY»NN _'J ' ] 1 

COMMON/ STOR 1 8/ TFR AC ♦ L I NE * PT 1 ME ,CONV » l ENP tlPR I NT » NPAG £ 

COMMON/S.TOR19/TIMF,TO’ ^11’ ' 

COMMON/STOR20/AREA, IOPT»TlM.'NTH»fHS»tHRUST 

COMMON/STOR23/SPH1 0 «CPHI 0 »LAMbAO >X O»YO»ZO' 

'COMMON/STO’R247'PHV,Mfn,NfAR»‘NPH»KPH ' 

C OMMON /ST0R25/L0W,FPTtNY»EPRIG*PJI,PI,PI2»DT,H I GH 

n T M r N s T ON _T T M ( 5 0 ) ♦ T hS ( 5 O) «S(9«6) *1^1 (.5 OJ » WPR ( 5 0 ) 

'P I M E N S I OJ _ NP AG E ( 1 0 L _ Z 1_ 

PIMPNSION PHT( 14rrNlp( 141_.MCH( 50J_tCPPi50JL 

D T MF^S I pN_S T ( 2J) ♦ 1 0 ) ♦ APRAY (10 *30 « 10 ) »SA( 10*2*10) 

real L at * long'* L AMiSaO * LOW 

REAL MCH, MACH _1 ] 

real ISM * INERT *MAS5[ 

external^ccflo' “ : 

PTTME s AINT('flME>Ptl ) »PTI -f PTI 

TLIM = PHT(21 - 1.0 

_T PR fLV - .= „T I M F 

CALL ATMSPH 

MACH = VT/SOMNO 

0 = (nENS*VT**2)/2* 

"call TABJ ] ^ 

PRJkG = cnS*o ^ ^ 

CALL "tab 2 ^ 

CALL TAB 3 _ _ •_ 

M ASS~'= we ight / GO ^ 7 ^ 


• OM»» » 


c ' tn 't f i al"i z'f la uk'ch varYablfs 

1 AXL s" (“thrust - *^[^PAr,) »TdT( fHFTA)/MASS 

A Zl" = ( THRU ST - HP Afi ) *ST N ( THFT A ) /MASS - (^0 

^ \jrHnn = 

_ _ ■ _ I F ( V t ' i GT’._ 1 . 0 F~ 5) GO TO 80 

_ ~ Dpf * VXL 'Z ~ 

VX L = O P T n«6D~A)*AVL 

DPI = VZL _ 

^VZL = DPI + Ti«6b- V)*AZL 

__ TIM^s TTME ~4'ud F-6 

“ ftO'CONTlMUF _ ^ 

VL = SQRT('VXL**2 + V7L* *2) 

c Tnttiali'zp S'a'pp'ay “ 





S(1,N) 


XL 


• 




S(?,N) 

S 

ZL 






S(^,N) 

s 

thfta 






S(4,N) 

s 

VXL 






S(5,N) 

s 

VZL 






S(6«N) 

s 

6.0 






S(7>N) 

’s 

AXL 






S(8,N) 

s 

_AZL 






.S(P*NJ. 

= 

THDD 






DO ,9P. : 

ST(j>i i « sjjvij : 

. oO CONTINUF 

lOO COMTINUE 

CALL PUK'GFl ACCFL.QJ 1 

T RDT .. =. T im ± S>1 

IF((TPDT PT.IMFJ.. .GT,..L0W._*AND. .TIMF.*LT.v PTTME) DT = 

1_PTJMP “ TI HF LO W : 

_ AXL_= Sl7*l) 

AZL.= S(8»1 _ 

VT_= SORT.(VXL**i +. yil*JL2J : 

CALL maTLCH 

_ XT = XL*CAZ ^ 

Y.t_= JLL* SJLZ 

ZT a - ZL • _ 

ZZT * ZT - RE 

CALL _TRAN_S2 ( XT„.YjjJ*ZT ♦U^.UX*DU; 




i 

i 


Y s UY 

2 = UZ 

A- =-_s 6 r II X * * 2_+ ~Y** 2 •<■ 2*» 

A LT = . 

VX, =.VXL*CAZ 

VY s VX L± SA2 

V Z = - V7L 

.CALL TRAMS?(VX»VY,VZ^h x,|Jy,o?) 

xnoT.,..=. -OMFr,A»v + i.ix 

Y no T _= n^"ECA»X + I IY 

2 HOT = ~ uzj 

AC X = A X LJ! CjA 2 ' 

_ACY =_AXL*SA_Z 

ACZ =._-AZL1 

.CALL. TRANS2 ( ACX » ACY » ACZ ♦ UX »UY , UZ ) 

._._XODOT_= UX - 2«*OMFGA»YnO f'+ ' 0 ‘mEGA;»*2*X_ 

YnooT = uY + 2 .*om> 6 a*xd 6 t +' omega** 2 *’y 

ZD Dof = UZ" 

CALL_0AT&1 

NN =_ 1 ^ 

^CALL sets A _ ^ ^ 

IF{ALDH *~Gt,“0.00 GO 10_75^ 

.TFRAC _'=...SA ( i7,2_» 1 JVjIa ( Uj *1) -'~SA(1 » 1 » 1 f) 

TFPa^ = fiPS(TFRAC) 

_ _’_TT'wF s TpREV"4^’Ypr&C*(TT*^F - TPRFVl 

XL = XL2 + TFRAC*'(XL - XL2) ' ' 

ZL'_= ZL 2 +”TFRAC*(ZL - ZLZ) 

_1_T yXL .* VXL2 + TFRAC*( vVl - VXL2 ) ‘ 

VZL = VZL2 + T F RAC* Yv ? L " VZ L 2 ~ )~ ' ' 

GAMA = a’ta'NzTvZL.VVL) 

EL = GAMA 

' GO TO 200 7 . '■ ■■ 

75 O CONTINUE ‘‘ ■ ■ 

!F(TImf .lT, PTT'm.fY gO T6“865 

^ imiME ".EQ.'Pf IMF) on TO'SOd 

C ' COMPUTE INTFPJ^OI.AT'tON FRACTION 

_ _ T FR AC”~="( PT 1 MF )'/ ( f I MF - TPRFV ) 

GO TO'sOf - ■ - - 

■ flOO tfrac's r.o 



1 



1 


PTT ME = TIME 



bOi 

CONTINUF 




c_ 

TALL PRINTOUT ROUTINE 



C ■* " 

w 


CALL OUT • ' . 



1 

c 

SET NFW PRINT TIME 

. 




PTIME = PTIME + PTT 




o 

cc 

i CONTINUE 




c 

STORE PRESENT VALUES OF ALL OUTPUT DATA 



NN = 2 

CALL SETSA 

TpREV ^ TIME 


XU = XL 

ZL2 = ZL 

VXL2_f_ VXL 

__ VZL? = VZL _ _ 

IFdlME ,G“T. 20*0) GO TO_200 
1 F ( T I w.E .* LT ,_T L I M) GO TO 100 

?oo return ' ' ~ ' 

* fnd'' 





IAtMSPH SUR, _ATMSPH ---^ATMOSPHPPJC VARTABl FS. 1 

_ __SUPROUTINE ATmSPH 

( “ COMMON/STORl/ALTtTFMP,PRFS»f>FNS«VTSC»FQUNn “ 

_'2_ TOMMON/STORfl/OMFGA ,RF,RA ♦PR,(;0 ,f,M,ALTM _ _ ' 

COMMON/S~TOP?6/PRFSO »ncNS0 tTF^^PO ♦SOUNDOVKl%COPr. KPP*PV 

( *__1 _COmmOn/STOR37/TSLOOF,mALTB,TFMP3,PPFSR 

^ pTMFNSTONJCSLPPFj _5_) tHALTBC _5J *TFMPR(_RJ_^PRFSR(_qj 

RFAL Kl__ 

( HALT = ALT*(Rf 7(RF + ALTn*rONl _ __ 

C ALL T_NjFRP(H^LT»N.5Oo<J,$80b)_ ^ 1 1 ~ 

DHALT = HALT'-”HALTB{M) 

( TEMP = TEMPB(N) + TSLOPE( N)*DHALT 

IF(TSLOPF(Nn400,500,400 

aOO continue 

( ^ PR ATlO = _PRFSR(N)*(TFMpB(N)/TEMP)**(KmSLOPF(Nj ) 


( 


GQ.T0.600 . 

(( 5PO_CONTINUE-. .. 

PR A T t O _ = ._PP.F$ B.( N ) * FX o ( ( - K 1 / T FMP R ( M ) ) * OH A I T ) . 

' A 00 continue 

( _TRATlO = TFMP/TFMPO 

OFNS_ ^0FNS0*oraTTO/TpATI0. 

SOUND = s”0UND0#S0RT(T9ATIoi ' ... 

C PRFS = PRATIO*PRESO 

_ _RFTURN 

r'OO PRFS = 0,0 

C DFNS' =~0.0 

DENS =.0.0 

^SOUNO = 0,0 

( RFTURN 

oPP PRFS = POFSO 

=' DFNSO ' * 

; _ _JEMP ’= TfmpO 

S OUND ~= SOUNDO _ J 

return' 

( : . ... . z. 


o 


V 


C 






t 


I 


( _ - - — 

&Cc}:kJ)ata _ _ — 

COMMbN/SfOR^e/PP'ESO ♦fS'fNSOVtFMPdVsoU rK'ppFv ~ 

( COMMON/STOR37/TSLOOF,MALTB»TFMPB*PRESP _ _ 

PJW^NSlONjrSLOPFLPJ , halts ( 9J .TFHPB( Q) ,ppp5P( q) 

^PfAL.Kl ...! I — — 

( data, TSLOPf / -6. 5F.-.3 »0. 0* 1 .OE-3f 2,8F-^ *0.^0’»-2'."bF^ ,-^.OF-3 »0 .V* ' 

1 0.0/ - 

DATA HAL.TB/. O.O.U ,bP+3,2P,0F+3»32*nF+3i47,0»i+?*52,0F+3,6r,7rF+3T' 

( 1 7®* E+3*8fl*743F+?/ _ _ 

OAT A J[f MP8/_288 ^Vs . 270‘.6A ♦ 270 .'65 .'282 .“65 . 

L__180_.65» 180,68/ " J 

( data press/ .1 ,0 , 2, 23361®-:! * *»,40328F-2 .8 .AAsS^F-a ♦ 1 ,0R488F-3 . 

l_l*A2.289f-4.1,797l8J-4.1,074JE-8»l,6223F-6/ ' ‘ 

DAT A PP E SO / 2 1 1 6 . 2 / , DFNSb/ ?'76RE - 3 /'. T FMPO / 2X8,18/. 

( JL SOUNDO/Xl 16 , 45/_. K 1 / 34 , 1 6 3 1 94E-3 / .CON 1 / , 3048 / , KPRFV 7 1 7 * " ' 

ENO 



1 

I 


i 

I 



h 


9T 


A>. 



cdatai sub, DATA1_---..C0MPUTP additional output 

3UBPOU.T I NE-U/aA] 

. DOUBLE PPECISION .CCS 

CO^^MON/STORl/AL T,TPWP, PRES f DENS., VI.SCt.SOUND. 

.COMMOM/«.TOR2/X,Y,7 ,R...J 

CnMVOv/STOR7/XDOT,YDOT»ZDbT,XPiDOT»''DDDT».ZDDOI 

COMMON/STOR8/OMFGA,RP,PA,RB,Gn,OM<ALLM 

C / S T ORJiZY L,y vL,'9XY,GAM L>OLDbT 

COMMON / S T OR n / X L » Z L » T w P TA» V X L , VZ L ♦ A L P H , G A M A , C A 2 ♦.$ A2 

_ _ COMMON/STOR15/CiVX,VY,VZ ,VT,EL»A2»RANGE»l AT , L.ONG, AC_CjGOLAt_ 

_ _ COMMON/STOR1<)/TIME,TO__ _ __ __ 

COMMOn/STOR 25/LOW»EPTTNY,EPBIG,PTI ,PI ,PI2,DT tHlGH 

COMMON/STOR27/I iPLAT,TlPLON,TlPR,TIPTlM _ 

^ _ .Cp_M_MpN / STpR.2 AilRLCH , A ZLCH «ELLC H, VfLLCH , F LCH « NL C H , 2 LCH , F D'l'CH , 

1 NblCHtZnlcH .... 

COMMON/.STOR 2P/.RPAD , AJ e AD_* FLR AD.» VFLVaO fFRAO » NR AO » ZR AD , E’DR A'Oi. 

1 . NDPAO,Zf^RAn ^ 

COMMON/STOR70/‘p)>LT*a »velVacc.t ] 

.real NRAr>»NLCH,NnRAn,MOLCH 

P EA L IIPLATt llPL ON , U.P P,I IPTIM 

.REAL lamda,mass»<,low,long»...lat.»..vach,mch»lamd.aj3 

VEL.E .SPpT(xp6t**..2 +_. YDOt** 2 _+.2bOT.**2) 

ACCi * SQRT(XDDOf#*i + VdDOT** 2 +. 2'6boT**2.l 

.„EX__=..XDDOT_ + .2,*0MFGA*Y.D0T__- .O.M;FrjA##2fX 

EY. .* . YDDO.T .2.,.*0MEGA*XDQT_ r._O.^EG’A*»2.*Y 

EZ^«.,ZDD0I 

. CALL, TRANSl(EX,FY,.eZ, ACX ,AC.y,AC2.) C. 

ACC . .* SOP.lf .{ ACX.**2 .ACY.*i2..J:_ACL*t2J 

CALL maTo.AD 1 

CALL TRANS1(X»Y,Z,WX,my»W7.) 

ZRAD ».-W2 . - : 

_ER A o = WY 

NR AD = WX __ 

RRAD = S6RT(VRAb**2 +lNRVbjL*2 j-' 2RAb**2 L I' 

_.CXDT * XOOT + 0MEGA*Y I 

CYDT * YDOT. rL_OMEGA_*X 

czbi • 2D0T._, 

C.ALL_r.R..A.N.S U.CXJlLt^lILiiiLQT.iJMX.iWX».yLZJ 

EDRAD * WY 

NDRAO,.». .WX ' ' I 

.. 2DPAO.».-WZ 


I 


-VELRAD„*„SOFJiFDR4DJL*2^m?AD*Jt2-+*ZD.*LAD«?2.1 

SRT » SOPT(fraO«* 2 + NRAr>**2i 

ELRAD.« ATAN2(2PAn,SRT) 

AZRaO = ATAN2.(ERAD.NRAD) ; 

IF(a£RAO .LT. 0.0j _ AZRAD-* .A2RAD..+...P12. 

-CALL MATlCH.._ 

CALL-IPvA^15LLXxL»L».h'JLil‘.'y.i.WiLl 

ELCH * -Wy 

-_NLCH-.».WX 

2LCH «.-W2 r;_RE. 

RLCH SORT ( ELCH*«2.. +.-MLCH**2-.+- ZLCH»*2I 

C ALL . TR AMS I .( aOT » .CXDXt C2 DT »>/X » WY > W2 J 

EDl.CH_^_W2T : 

ZOLCH s -WZ ■ ■ ■ 

NOLCH.* WX 

VELLCH * ,SORT(.EDLCH**Z +„MDLCH**2_ + ,.ZDLCH**21_. 

. SRT * SOPT(BLCh)**2 _+_MLCH.**2) 

- . .ELLCH * ATAN2.(ZLCH,5PJ1) 

A z L CH ..A IA.N,2i£i_CH J.N.LOJJ 

IF(aZLCH .LT, Q,0)_AZLCH__« AZL.C>!_+_PI2. 

XL =-ELCH*SAZ„+ NLCH.*CA2 

YL .»-ELCH*CA2_.+_MLCH*SAZ 

JL- = ZLCH 

yxL Eo.LCH*sAz. . molch^caz 

IV YL ^^EDLCHtCAZ MOLCH^s l Z ^ 

VZL « ZOLCH 

RXY = .SORT(XL**2 + YL«*2J- 

... . SRT * SQPT(VXL**2 +..VY.L*«2) 

6AML. * ATAN2(VZL»Sp.TJ| 

RLOOT s (XL*VXL -+-YLtVYL..+ .ZL«VZLi/RLCH 

^pAKGE_»_Oj( 

1F(RLCH ,LT,..5.0) fiO-TO.7.00 

. . CCS * (R**2 + PE**? - RLCH**2)/{2,*P*REJ 

IF(DABS(CCS) .GT. 1 .OoO.) CCS * 1,000 

0ELTA..« DATAN2(D$0P.t(l,-CCS**2) .CCSl 

RANftE * RF*OFLTA.. 

.2 0 0_CPMJJ.N.U.F 

£RT » SQPT(Y**2 +Jt**?) 

PHI ? ATAN(Z/SRTf . -"I 

APfr J*_j RA/RR) *_*21< zZSPii 


1 


! 


f 


V 



GOL&T = ATAfi(ARill ! 

LAT = PHI ; 

.... LAMnA_ = . aTAN2(Y,X) 

LONf, i LAMDA __0V.cr,A-»(TlMF - TO) 

IFaONG .GT, PD^L^Nr, j= LONG - PI2 1 

L_ I FJ UOMG • LT „(_-o.T )) = LONG ~ + ^012 

FL = ATAM(- VZ/S0*?T(Vy»»2 + VY**2 ) ) 

ALPH = t heta - EL 

. A i . = A T A N i.( ,V Y ♦_V X ) 

.J.F1AZ . .UT . 0_. L.a2_=_A2-J^_P.I 2 

IF (ALT ♦LT, 5.0 . AMO. ^L .LT. 0.0) RFTURN 

.._.CALL UP 

RFIL'M 

FM) :l 



I 


I 


’ % .JL 


T 



COIRECT SUB,. DIRFCT T HRUST CONTROL OPTION 

__ _ subpoutine direct _ _ ' - 

COMWON/STOR 1 /ALT ♦ TEMP , PRES /PENS »V I SC »s6un5 

COMMON/STORn/THCOM»Tl,ALTC*THX,THY»THZ 

COMMON/STOR13/XL»ZL»ThETA,VXL»VZL»ALPH*GAMA,CA2,SAZ 

.. C0MM0N/ST0R15/Q»VV.,VY,VZ >VT.»EL»AZj».RANG.E»1 AT,LONG»acC*GDLAI^ 

C CmmON / S T OR 1 9 / T I M F , T 0 

C OMMON / STOR20/APPA, TQpT«TTM « NTH*THS» T HRUS T 

_ dimension TIM(5Q).,tHS130J 

IMTEGER.THCON 

EL . ATAN..{~V.2Z.SQPXLVX«A2_+_.VY±*^^^^ 

XFRAC .?...VX_/VT 

yfrac_*..vy/vt 

2F«LAiL-= VZZYT 

CALL. TRANS2 ( XFRAC ♦YFRaC»ZFRAr«UX» UY« U2 ) 

XFRAC. =.UX 

YFRAC_=_.U.Y 

ZFRAC “ _U.Z ^ 

-THFT.A =_E.L. 

IFdlME .LT . T1 ,OR« mSET «EQ« 1 «O R. ALT ,LT, ALTC) GO TO IQQ 

. XFPACi = XFRAC 

YFRACl = .YFRAC 

ZFRACl = ZFRAC 

..ELST .= EU,. 

TIMST_= TJMF 

.N s E.I.. = i : - 

.jOO CONTINUE ^ 

. IF(TIME...LT....T1., .OR..ALJ.. •.LT.„ALIC) ..60_T0_. 2.01..' 

60 TO {2Ql.»202)^TWCQi} 

_.?Cl CONTINUE 

C NO THDUST .CONTP.OL 

THx - E . _T hr us T *X FF Ai:^ 

. ..thy = THRUST»YFR.AC 

TH2 = THPUST*Z£RAC 

60. TO 900 

20? .CONTINUE 

C-... THRUST CONTROL 0Pf»0~N . 

_1 'I’HX = THRUST»XFRA C1 ' ~ ~ " ' 

t*HY = THRUST^vYra'CI __ _ 

" THZ = THPUST*ZFRAC1 ' . 7 V..'.-” 1 . 

call TRANS1(XF.pac1,YFPAC1»ZFRAC1»WX*WY»WZ) 

. THFTA =.ATAN.{-WZ/S0RT(WX»*2.+„KYX^2.).I 

oOO RETURN • .’ 



Cl TP _ <5^ B •_ I I ANTA Nf OUS,J MPACT. POj NTS 

SUBPf^J I NE_U^ ~~ 

_ dou'rlp : Decision s^iri?‘,csiG2.siG»SQ _ 

DOUBLE PRECISIOM. FCE » '^£‘^A?CE I P , FSB I P,» SDEL ♦ CObL vDEL.*F.».G*A 

DOUBLE PPECISJON ...X I P »X1P»Z JA.»LjC_» ESO ♦XO^* YO » ZO »X ♦ Y,»2 

COMMON/STOR2/_ XX ,YY»ZZjPP. 

C OMMON / ST OR 7 / X DOT. vYDOj , Z DOT , X DP.CT .♦ YDDO T » Z DDOT 

C OMMON / ST OR 8_/ 0.^1 E G A »PR*r- 0 , G Y , A L I M 

C0MY0N/ST0Ri3/XL»ZL .THET A ^VXl ,'v.2L>AL^ »GAN'A.,CAZ»SAZ_ 

COf/MON/STORlR/TIvPjTQ 

COMMON/STGR23/SPHIO»CdhI b»LAYDAO»Xx6»YYO,ZZO : 

C OMMON / S T OR 2 5 / LC W , F P T T N Y » cP B I G , ? F I , P I , P I 2 , tn » HIGH 

_ COMVON/STOR27/I IPLAT, TJPLON*! IA>P»I IPTIM . 

^AL_L Af'^'l^lAS^ , K, LOW ^ LONG ♦ LAT, MACH , NCH > LA VPAO ,LBAR 

real. 1 1 PL at , r I PLON ,T I_r>R ♦J.J P.T I M ; ] ■ ■ 

C I N I TJ A LI ZE -_.X , Y ^Z..* XDOT.-. Y,DOT ♦ ZDOJj 

X_ = , XX 1 

Y_ =_YY : _ 

_ _ _z = zz : 

R_-_RR_ 

’ ySO =. x'p6t**2 -I- Y D0T*«2 + ZD0T**2 ~ j ^ 

SRGM .= SORf(GY)__.__J 

ECE = R*VSb'/GM -_1.0 

A = R/{1,0 - Ecn _ _ __ • 

ESE _= { X*XpOT_ + Y*YDp_t__+_ Z*Z ppJJ / ( SRpy^pSORT { A )J 

FS O = _FCEAf 2__+_R§f*t2 ^ 

E c F I P = . 0 OO _ P F z A ; _ 

.SO = ESO.- fcfip*«2 : 

_IF(SQ .GT. O.OdOT go to' job 

so = O.OnO _ ” * ^ 7 J 

100 ESEIP = -DS0RT{s6)_ 7- ' 

SDEL. B. .'/.GEJ P*LCE ,^_Frf I_P*E5 E ) /E_SQ_ 

CDfl = M^CEIP^FrP + pS5J.P»‘",SF.) /.FSO : 

._..DEL = . DATAN2{SDFL,rD,PL) 

IF(DEL .LT. O.ODO) DEL = DEL +_PI2 

F = (CDEL - fcp)/(t,o -_«^CF) _ 

G = (A**l .5/SRG^M^^ (S<^PL + fsf « FSFiP) 

' Tf'l = ( A**1.5/SPG^)*(nE 'L V ESF FSEIP)" ' _ 'I__ 

CA'.-'T* = TEL + (Y IME' -"Tb) “ “ ■ 

XIP = F*X + G*XOOT ^7 ” _7_^ " ' r Z7_ 1 

YIP = F*Y + G»YOOT 



j: jj? _=_F»z dot : 

LBAR .= LAMOAO + OMCGA»CA PT ■ 

XO s.RE*C°HlO*COS(LBAe) [ ” 

YO = .RE*CPH!0*SIM(LBAR) 

ZQ r_RF*SPHlO _* 

LC. =PSQRT( (X!P - X0)»*2 (Y!P 4 (ZIP -Z0) »* 2) 

L A IA.M2 1 YJ P.,Xlij I 

LONG ..= ,.L AMDA. - ,.0'^EGA*CA.P_T 

PHT =-ATaN(ZI=>/50RT(_XTP»*2_+^YJP**2) 1 

GHJ^Hl =. . ATAM( (RAARP.l*»2*ilR(PHn/-COS(DKI U 

SSTG2-= LC/(2,*REL 

I F.( dabs ( SS I G2 ). . GT l.ODO L^SJ G2. I..QDJ2- 

CSlIl2_=„.DSGR.Ij.l .0„^_£SIi5L2.»»2J 

SI6 = 2.*DATAN(S_S.1G2/CS_LG2J. . 

RANGE =.RE*$.I5 

II PLAT ^GDPML 

JI^LON = LONG. 

II.P.R = RANG F 

U PT I M = rAPT ■ 

RETURN 

.E.NQ I 



NiPUT SU9, Input rbao input ano_printout 
subroutine input _ „ 

dimension AA<14) J 

DATA EN/6HZ99999/ 1 

rewind 5 _ _ _ 

23 continue r_ _ _1 

READ(5i_19)_AA ^ 1' L_ ___ 

l'9 F0RMATC14A6) 

. PRINT l9, AA 

IF(AA{2) ,EQ.. EN)_Q0_TQ 34 

GO. TO. 23 __ _ 

34 REWIND _5 _ 

.return 

end 



I 


( 


_ CTMJPOP. _ . _ 

( ’ SUBROUTTMF TNTFPP(OG,>fP,*,*) 

~ "cOMMON/STOP^6/PPF50,ncNS0»TFMP0»5OUNn0*tg»CO^1 »kPRFV 

C0M'M0N/ST0R37/fSL0PF.HALTB*TFMPB.PPFSR 

( ' _ ^DIMFNSION TSL0PE( o)»HALTB( P),TFMPB( 0),DRFSR( O) 

_ KP_= KPRFV _ „ _ _ 

iF(OQ .LE.‘ HALTB(KP+1) .AND.* 00 .GT. HALtB(kP)) GO TO 2o6~“ 

( ■ * “IFC'OQ .LE. HALTB(D) RFTURNl ‘ ■ " _ ' _ 

"If ( 00 .G'=’. HALTB(9n PFTURN2 ' * 

IF(0Q hprfv) T30,20bVioo 

( foo j)0 jl ■ j=KP,8 

IFIoO .GT. HALTB(J> .aNH. do .LF. HALTBt J+1 ) go to 120 * 

ilO CONTINUE ” * 

( 120 KP ■ = j * ■ * ■ 

I * GO* TO 200 ' ' 

l50nCPP = KP - 1' ' " _ 

; ( - - - do“i6 J = 1»Kdp ' ' * 

: ■ "iCPJ = KP - J ■ - - 

_ JiF'lOO .GT. HALTBUCPJ V .AND'. 00 .LF. HALTB(KOj+l )') GO' TO l 

C { __ 7 J60_CONTINUE - - 

! ■ 180 KP '=* KP - J 

' f 7 00 CONTINUE ■ 

; C -■ KPRFV = KP 

' 'H PRFV = 00 ^ ' 

__ RETURN ’ _ _ * _ J"' 

• ( 1 _ _ _ 


( 


C 


( 

i — 

' . c i: 
? 


( 

• • i 

} 



CO 



C.LaUMCH . . Slip, LAUNCH LAUNCH RAIL . • • _ 

surpoutime launch ] 

DOUBLE PRECISION ST . 

COMMON/STORl/ALT,TPVP,PRES»nFNS,VISC*EOUNO ^ [ 1. 

C0MM0N/ST0R4/SARPA,MrM.C0n*O.»MACH»NCn»CDS»0RAG 

C0HM0N/ST0R5/TWT«NWT,'.-'PR»WEIGHT,MASS . _■ 

''_rOMMON/STOR8/CHEGA,RC’,PA ,Re*Gp*_GM_»ALfH2'2" __1_1 I 

C O^^MON /J, T_OR 9> ^ MK - 

CO^^MON/STOR'iO/ST _ _ . 

COMMON/STOR15/Q .VxTvV ,VZ »VT ♦ eX *AZ *RANyf»l AT ♦LONG, ACC *g6j[^^ 

COMVON/STORl7/LENGTH,S E L_,_D I_S T ,_L^ L_, L A C C_ ____ 

_COM'^ON/STOR1R/TIM|?,TO 

COMHON/STOR20/APEA,ICP'f ♦‘TIM,NTH‘,fHS,tHRUST 

C OMMON / SJO_R 2 5 /L0W,EPTTNY ,'EPB IG,PTI ,PI ,>12 ♦D T'/HiGH ] 

DTMRNsiOM' ST(>0,10j ^ ■ 

dimension MCH(.50) ♦CDOCSd') ♦ T WT ( 50 ) ♦ WPR ( 50)* T I M ( 50 ) ♦ t~HS ( 5Q ) ♦ S ( 9 ♦ 6 ) 

real LVFL,LACCXL'^\X,jH,LviL2.,LACC2,MASS ”_1 

FXTc-Rmal.ACCEJL? 1 

IF ( Lf^NGTH ,LT« d.ll ocTL IRN 

_C A LLALM S P.ij I 

TIMP = TO ^ j 

X OT = AMiNliP,p.l_,.0JJ 

CALL .T.AB2 

CALL..TAB_A : ^ 

MASS .=. WPIGHT/.GO 

L AC C = (TH RUS T - WF IGhT»SFL1 /MASS 

I ' s ( r, N r7.~Q7'ii ^ XL ■ ' ■ ■■ 

S(4,N) = VT 

S(7,N) =_LACC , 

ST (1,1) s S(l,l) ' 

_ST(4,1 )__=_S(4,jj 

ST ( 7 ♦ 1 ') s ST7 Vi ) ^ 

1 00 CONTINUE 

D I Sf 2 = p 1ST ^ 

ALT? * alt 7 

LVFL2 = LVEL 

TiMf:2 *_TIMF_ ^ 

CALL PUNGF(ACCFL2) ‘ 

DT = AMINI (0.01,DT) 1' _ ' 

ALT = OTct*SPL _J__ 

...IFIOIST. .LT, LENGTH) r,o TO 100 

TFPaC s (LENGTH - ^IST2) / (DIST - OlSf?) 

ALT = ALT2 + TFRAC»( ALT, -_ALT2 ) 

^ — LVEL = LVEL2 + TFRAC*(LVEL - LVEL2) 

.. time = TIME2+ TFRAC*(TIME-TIME2) 

VT = LVEL - - — 

„oOO_.RETURN 

END 





.. -CMAT ... . .SUR. MAT_ --- ROTATION, MAIPJX. ANO INVFRSE 

.....subroutine .mat : 

DOUBLE PPEC.ISION, TMAT»SWAT 

. DIMENSION. TMAT( 3»3 ) ,..RMAT.(.3_, 3) 

COMMON/STCR.3/TMAT , .SMAT>S.PJHI.».CP.HI.»CL»SlI 

IMAT.( Lj .1 -L- !L_-S PHJ *jlL 

TM ATI 1 ♦ 2.L5 _“E PH I *%L 

that ( 1 ♦ 3 ) . .= . C.PHI 

TMAT(2vU_=--SL : 

that ( 2 ♦ 2 )_=_CL 

. TMat( 2»3.)_ =...0. 

I M A.Ii3_».l ) = -C PHI»C L 

that ( 3 »2) ..= . -CPH.I.’LSL 

TMA^ ( 3 ♦.3J_..=_.-SPHI 

SMAT(1.11_= TMAT(1,U 

SM A T 1 .1 ♦ 2 )_..= ...TMAT (,2,.1J ^ 

.SMAT 1 1 J.3L. =_J.MATJ.3 J.U 

.5MAIJL2j l.L_=_IflAUI.,2J 

SMAT(2»2) =..JMAT(2.2) 


SMAT(2.3u=-TMAT(J3,2.) 

SMAT (3.1) _= .TMAT(1,3.) 

SMAT ( 3 .2 )_.=-JMJSi.T(-2.»3J 

SMAT ( 3 .3 )._ .r-.TMAL(.3 ..3.1 

RE.TUJR.N 

END..._ 



CM&TLCH .. SUB. MATLCH_--r .LAUNCH-SITE_MATP1X_ 

SURRQUIINE-^IAILCH 

DOUBLF PRECISION TVA.T.SWAT 

common/stor3/tn!at,$n‘at»sdhi .cphi .cl* si 

COMMON/STOR8/OMEGA,RE*FA*RB.GO»GM*ALIM 

COMMON/STOR19/TIME*TO_ . 

, C OMMON / S T OR ? 3 / S P H 1 0 > c p H I 6.* L A M DA 0* X 0 * Y 0. * 2 0_ 

D liil E N S 1 AT(3*3 WSMAK3.31 

REAL._LBAR»...LAf^DACj 

SPHI. -..SOHIO 

-_.CPH1 .? .CPHI.O. 

LBAP = LAMDAO + OMfga»(T!M E - TO) 

-I C U_=_XO.$jiMEU 

CALL MAT 

RETURN 

END ^ 





0>>w c 4A& 


TR AH ■ s"u'r . ” ^&J.R a'o 1-~_ 'ra HAR ' SI fp~'MA tRW 

S UR R 0 UT I M M A T^A P) " " ' 

DnURL«^ PRFCTMOw' TMAT*SMAT 

_ . ” . common / S T or R / T m A T , S m A j ,, 5 p'h fvC PH T » CL » Sj. 

C OMVON/STOR8/OMEGA.PE , PA» .? Bj.G.6 i j. A L J M__ 

COmmON/STOR19/TIME»TO 

COMMON/STOR31 /SPHI 1 ,CDHI 1 »LAMDA1 

OI MFNSIONI SM AT(3*3^ *TMAT(~3«3) 

RFAL. LRAR ,LAMD.A1 

SPHT s.SPH.U 1 1 

CPHL =„CPHI1 ^ ^ 

L 5 A.R „ = ..LA Mp, A 1 + O VFGA »(TIM g ~ Ip_l 

SL..r.._SIN(.LBARJ 

CL = C OS (LBAR) 

CA.LL. MAL I 

RETURN 

END 



CMaTROC sub. L.PPFSPNT ROS.ITIOW MATRIX 

..SURROUT I ME-MAJJ?Q.C 

.OOUfiLP POFCTSION TyAT,SMAT__: 

C0MM0N/ST0P2/X»Y.Z*R ... 

COMMON / ST OR 3 / T*' A T ♦ S V A T ♦ S PH I ,.C.PH I ».CL «SL 

_ OTM.PNSION TMAT(3.»3 L»_SyAT(_3»3.) 

R ..= . SORT ( X**2..t ; 

S RT = SQ PT(X*»2 + Y*»2~1 

SPHI.-= .Z/.R 

CPHI =_SRT/R 

_..CL..» -X/SRT__ 

SL..= .Y./.SRI 

CONT.INUE 

CAU_„*4AI 

.R ETUR N L 



I 


C_QUT. 1 ' .SU3. OUT INT'FPPOLaTe DATA AND PRINT 

SUBPOUT l ML OUT J ~~ 

_ COVMON/STOR16/SA »APRAy,NN _ _ _ 

COVMON/STCP 1 8/TFPAC tLINF^PT lMEt CON V » I F ND > 1 PR I NT » NP AGF » KP AG' 

C OMMON /.S T OP 2 1 / w F S S j PH A 4 

COM»*ON/STOP32/NSYS »COM ; 

COMMON/STOP33/''! i P,NCy.C.».NLIN 

0 OMMO N / Sjrpp 4 / < 1. 1 w 

6 1 y P N SION COn'( 1 O" ♦ 7 ) .♦ A PRFY(10/30 »10? 

_DlyENSlON MF5S(30) 

DIMENSION NPACiEdOl ‘ U„ L 

D I MENS I ON _S A.(.T_0 ♦ 2 » 1 0 ) «ARRAY n 0 »30 . 10 ) »N VAR_oi^l»__li MA ( 3Q j 

__REAL MESS , r 

^iWoPMAu iMi ) : 

,2 0 .FORM AT ( 1 X » A6 , 1 X ♦ F 8.2»«^X*F6«3 *5X*F8. 3»5 X_» F8 • 3 « 5X » F 8 • 3 ♦ 5 X . F 8 • 0.» 5X.» 

J Fn«2»?X»FE,4f^X»P0.4) 1“ 

_2 1 _ F pP ‘M T ( 1 X ♦ A 6 ♦ 1 X ♦ F 8 , 2 , F X » p_8_. 0 ♦ 5 X_,.F7 ..p.* 3J^»_F 7.* 0_» 5 X, F 6.* 2 ♦ 5 X . F 9,2 » 5 X * 

1. FP,2,EX,F10,?,FX»F_7.1 ) ... „ 

22 FORMAT(lX»A6»lX»Fe,2t7X,F9.2»2X»F8.3,?X*F8,3,?X,Fl0.1_*2X* 

i FQ.2.2'x%F9.2\2X»FQ.2,2Y«F10.1«2X»Fl0.1*2X*F'i^«l) 

23 format ( ix»A6»ix»F8,2.?X»Fl6.6»2X»Fio.6»2X»Fio*0»2X»F9,l»2X»F9.1» __ 

i 2X,F?,l»2X«Plb.0»2X»F7,3»2X»F9,l»2X»F9,i ) 

24 FOPMAT(8X»8H TI_ME _* 2 _X* 10 H XL .2X»l0H _YL _ »2'x‘, 

i iOH __'ZL »2X‘»0H_ VXL _ »2Xt9H VYl ,2X»9H VZL __#2X» 

~ " 2 IOh” _iXY ♦2X»7H 3 amC’ *2Xi' 9H'_ VEC-L " t'2X*9H RLDOT j. ' 

25 format ( 8X » 8 H ~ ('S F O » ?X*10H _(FT.) »2X»10H (FT) _j2X*^ 

'l i‘0h_ (FT)~ ,2X»oH (FT /SEcT» 2X.9H ( FT/SEC ) »2X ,9H~ (FT/SEC) f2X »_ 

_ "2 lOH (FT ) ' ,2X,7H (f>FG) »2X»9H ( FT/SEC) »2X»9H (FT/SEO_//) 

'26 F0PMAT(lx’»A6»lX»F8,2*RX*Fn.4,5X*F9,4,ex»Fln.2»5X»F8*l) 

27 FORMAT(lX»A6»lX.F8,2»ex»F11.0»5X»F8,O,5X»F8.1) 

_ _ 28 FOPMAT( 20Xt47HINSTAMTAMFOUS IMPACT POINTS (VACUUM TRAJECTORY)//) 

30_FORMAT(8X»8H JIME ^»5X»8H AL PHA »5X ,8H TH FL »5X»8H FL EL t 

.1 5X,8H FL AZ ♦‘^X^8H ^ALT...»5X t UH ..RANGE .. *5X»8H LAT._<3D f 

2 5X»9h long. ). 

31 F0PMAT(8X»3H (SEC) »«X»8H_ JOEG) fSX.SH (TtFG) f5X,8H (0E6) » 

1 5Xt8H (DEG) »5X,8JH LPJL_»5X_».li»£. 1N.M.)_ »5X_»8H _iOEG)_* 

2''5X.9h “" ('d’eG) //i’ 

32 FOPMAnSX.BH TIME t«^X»8H THRU ST » 5X»7H DRAG * 5 

"*T 5xi6H MACH ♦Ex T^H DYN PR »5x78H rTl’ ACC * 5X . *9H "REL VEL ']»5X'* 

_ 2 7H mass ) ^ ^ . i 

33_FOPmAT(8X»8H __(SEC) .sx.8H'_ (LBS)_»5X*7H (LBS)*5X,7H (LBS)* 







I. 


i'YiX,5X»PH(LSf7FT2 W^x“; 9H(FT/SEC2) »‘5X,10H' (FT/SECr »‘5X»“‘ ' 

_ 2 7H (.<LUrj)//) ^ 

_ 34 FORyAT<8Xi8H '* TIM'P , 2 x" H~’ S L' RANGE T2X^ 8H' LOOK J^ZT2X teH LOOX 

i_2X.10H _VEL ■ V2 ^,oh east “*2xV9H NORTH " ♦*2X»9H" VERT T 

2 2X,10H "vEL-_E' .2'X»10H ’ ^V^L-*N .2X,‘10H " VEL-V ) _ 

'?5 FORMAT{8X»8H ' ( SPC) » ?X . 9H_ (N.M.' V»‘?X t PH "'VnPGV *,2X »'8H‘ (r>EG) ~« 

_1_2X»10H (FT/SECJ ._2X»<5h" _ (N.H,)«2X.9H (N.M.) ,2X.9H " ( N.M. 

__ 2 2X,10h (F't/$FC) ‘♦9X,10H fFT/SEO ♦ 2X » llDH'^t Ft7TECT“/71 

36 FORmAT(20xV3]H launch "SiTF rOORniNATE SYSTEM*//) 

'■■37 FOR^'AT(20X♦>OH RAOar VitE' COORDINATE <YSTFM//)- 
38 FORMAK 20X»2PH I N^Pifl AL COOPDINATP SYSTEM '■//)'"■ 

42' FORMATi8x78H "tTmf" ,‘sxV'8 HLAf G6 »5X,9'H "LONG iSxTfo'H RANGE" 
•Sr i.8X,8H io IIME)" ■ *“ ■■* ■■■ . * ■ ■* ■■■■ ■■ ■ 

43*T0RMAT(TxV8H ' (S^C) »«X» 8H (DEGT »5X,9H fi>F<5T“»‘5 X » YoH I'Ni'M.) 

i»SX,8H'(SEC) _//) ^ _ 

a4 FORMAT ( 8Xf8 H T J ME__ »_AX , 1 IH R *5iY£8hJ' VEL_ »5X_» 

1 8H ACCEL ) ■■■ * 

45 F0RMAT(8X,8H TsEC) ,SX.11H (FT) rSxVeHtFT/'SEC) f5X»” 

_i oh(Ft/sfc2)7/) ' ■ * - ' - - 

_ 5 1_ format r8X»8*H ("YEC ) _♦ SX^SH _<DEG') .JX^RH ( DJG") #5Xt8H rOE'(jTT 

1 5x',8h ' (DEG) 'V5XiVh (M')~ ;5X\11H (KM) t'5xV8H" (DEG')","" '*' 

2 5X,PH .JDFG) //) ■ ~ ■■ ■_ "Z IT 

. --53_fORMAT(8'XjRH ( SEC )_ » *'8H *NT)' *5X,'7H' (Tif )’» 5X » 

__ . 1. llX»5X»9H_iNj/M2)__.4X,^" 1M/SEC2 ) . 5'x »'lOH^iM/SEC V 5Xi„„_ 

2 7M (KG ) / /) 

__55 F6 rmaT(‘ 8X»8 H TsEC) »2X,9H (KM) f2X.8H (DEG) ,2X.8H (OEG) 

1 2X»10H. (M/SEC) »?X, PH (KM) »2Xf9H (KM) »2X|9H (<M) f ' 

. 2 2X.10H (M/SEC) .2X,'OH (M/SEC) ,2X,lOH (M/SEC) //) 

_..56 FORVAT(8X»8H .(SEC) .pX.lOH (M) ,2XflOH (M) ,2X» 

1 10m (M) _ ♦2X»PH_ (M/5ECjf2X»9H ( M/SEC )» 2X 1 9H ( M/SEC ) 1 2X » ' 

_ 2 lOH (M)__ ♦2X»71^(ne_ro i2X,9H (M/SEC )» 2X *9H |M/S_EC)i^ ^ 

_P3 FORmaT( 8X.8H_ fS'^C)' fe^»8 H~JOFG') TYx",9H __^(DFG*) »5X.10H" (KM) 

1»5X,8h (FECi_7/i J - ^ 

05 FORMAT (8X.8H (SEC) '»AX'*nH_ JVm ) ' ' »>X ♦ 8)H ( M/SEC) ♦ 5 X'»“ _ 

1_ OH (M/SEC2)//) ' ] _* 

lOl FCPMAT('"i20x»6H~PAr,E ,.aIihA //1 ^ 1 II 

J02 format ( 120X,6H PAG^* , 14 » 1HR/>J I !__III___ 

..1O3 format'? 120 x'» 6 H oaGE VlAilHC//) _ 

104 FORMAT ( 120X.6H RAGE « 14 » IHO/"/ )' " ' ' 

_105 FORMATL 120X,6H PA_GE_Y 14 » IHF/"/ I Il_~ ' 




lO 6_ F0» V AT ( .1 ?0X , 6H =» A r,E J 4 «1HP/ /) 

DATA BLNf /0^77'^77777772/ 

data PLK</0202020202P70/ 


DATA-NVA7/8»8»10,1.0»4,3-' 

I F ( N 1 1 P . . EO . 0 ) „GO TO 3 0C__ ._ 

L . « .LINE_+J 



.. NUN * -NLIN ± 1 

IF(MLIN .EO. NCYC) Gp_10_j0_0 
GO TO 301 . 


3 OO 


F ( L I N E_. • EO . 30) LTN F=0 


361 


NLIN » 0 

-UXNF^,CjNL,t.l 

L_« J-TNE 

.iaiN_»_L 

MFSS(.LL.?.PHAS_ 

PHAS .= bLkK 

.CONTINUE. 



TTMA( L) « P TINE 

DO 11 J = 1,6 

NV « NVAR(J) 

06 10 I » 1,NV _ 

ARRAY(I»LtJ) » SA~(i, 2 ,j) ♦ TFRAC* ( S A ( I , 1 , J) ■ SA(i, 2 ,Jn 
APR EY t i L * J ) _ s_ ARP A Y ( Y VC , J ) *CON ( I » J } __1~_ 1 _ 

_ I F ( S A(.4 , 2 _j 2J_. GE ^ 0 Yo . AND,~ SA ( ffl', 2 Y' , G E > "o , 0 ) G0'~t6' 10 

APR'aY"(4~»L»2) » BLNF 2 ; 

AOocy( 4,L*2J » eLNY~ ~ ^ 

10 CONTImUE ^ _ 

11 CONTINUE * ' _ 1. T! " 

convert to OEGPFES 

ARR AY (.U L ♦ 1 ) = ARR Yf JJ ,L ,_1 I^CONV 

. _ ARRAV(2»L*i.l » ARRAY( 7,L,1 )/CONV._. 

ARPAY(3,L»1) - APPAY(3,L,1 i/CONV. 

ARRAY(4»L*1) * ARR AY ( 4 ,L , 1 ) /CONV 

ARRAY(7,L,1)'» ARRAY(7,L,1 »/CONV' J _ 

■ APRAY( 8 ,L»H' * ARRAYtR,L»l)/CONv“' Y' ' 

ARR.AJT (_8 » L ill g A RR AY ( P ,L . 3J /CONV " 

ARRAY(2»L*A) » rPRAYTtVLTAt/CONV 

ARPAY<3»L»4) = ARRAY /CONV ’ Y. 

_APRAY( 1,L»5) g_ARPAY(ljL»5)/C0NV7 ’ _ „ II 



w 


f 


( 


VXL» VZL - are velocity components in the downrange and vertical 
directions. Found by integration of AXL and AZL. 


CS 

THETA 


FNORM 

CNA 

SAREA 


COSfTHETA) 

is the inclination angle of the rocket, the an^e between the 
principal body axis and the hciizontal. It is originally set 
to the launch elevation angle, and is modified by double 
integration of THDD, THETA douUe dot. 

CNA«SAREA,i.Q *ALPH 

is the slope of the coefficient of the normal force acting 
on the center of pressure. 

is the reference area. 


Q - is defined above. 

ALPH - THETA - GAMA, and is called the **ang^e of attack". The 
program switches from the Period I model to the model 
used for the rest of the flight when ALPH becomes zero 
(found by interpolation as it crosses from plus to minus). 
The duration of Period I is about 0. 8 second for a 
NIKE-CAJUN, and about 5 seconds for a SCOUT. 

GAMA « ATAN (VZL/VXL) 


GAMA is the flight path angle which defines the dii . 
in which the rocket is moving 


MASS 

AZL 

THDD 

LSM 

LCP 

LOG 

INERT 


is the weight of the rocket divided by the force of gx^ lty. 

is acceleration in the Z direction, which is vertical at 
the launch site. 

is THETA double dot, the acceleration in the rocket's 
inclination angle to the axis. 

LCP-LCG. 

is the distance from the reference position (usually the 
nose of the rocket) to the center of pressure. 

is the distance from the reference position to the rocket's 
center of gravity. 

is the rocket's pitch moment of inertia. 


The equations above have been simplified bF the elimination of effects such 
as pitch damping and jet damping which were found to have no significant 


r 



ARPAY(2»L»5) = ARRAY(?>L*5)/C0NV 

ARRFYYI » Li l'r~= ARPFY ( V,L » 1 )‘/CONV 

ARRFY(2»L.1) = APRPY(2.L»1)/C0NV ! ~ ' 

_ ARPPY(3»L»1) = ARPFYC^,L»l. )/CONV 

“ ARRFY(4»L*1 )’ =' ARR^Y ( A ,L ♦ 1 )/CONV ] ~~ 

^ARPPY(7.L»l)"'y ARRPYr»,L.l)/CONV_j 

ARRPY(8-«L»1) = APRPY( A»L »1 )/CONV ^ ^ 

ARP cY ( fl » L ♦ 3Y~= ~ a PR ^Y ("P *L"» P ) >f ONV 

^ ' ARPFY(2»L»4) = ARRCY( ?,L»A)/CONV _J 

_1 " ARR‘^Y(3.L»4) = APRFYV^'.L »A) /CONV _21 ' _ 

Il_ __ ♦L_»F ) = ADPcy( 1 ,L»M/COKiy ^ 

l’_ FYJ 2 ♦ L ♦ 5j/ Y'A RRJ-;5^( Y /CO NV ; ^ 

_J*i____CAYL _TP01JT • ^ 

C._ i"F LiYf = 3'0 op IcNn = l OR I PR IJl’’’’ =_i » p|!intout__ 

_ J _ '_‘l F ( L I N E , EQ'.' '30 Y'O P r T FNO VFQ 1 « OR . j £R 1 NT . EQ« 1) 

GO TO 200 J '3 _ 

2_’*"50* CONTINUE ^ " | 

' kpag =_VpYg + i" _ ’ 

GO TO (14l if52,r51 ) ,NFYS ^ 

__151 CbNf‘lNUE'___ 

I . =-''i_ rZLZZ I ■’ 

NPAG_=-NpAGFLIJ 

I =_.I„ + .1 ^ 

GO TO (61 ».62,63..»64»65.»A6.»Y7,68).»NPA5 

6 1_..C Q.NT J. N.U E I 

...PRINT 101.,KPAG 

print. 30 

print. 31, _ 

C . .PRINT page a (cmglTSH) _ _ ... ... ..... 

PRINT 20,..( (MESS(LL) ,T.IMA ( LL {.ARRA.y’(.I..»L.I_,.1 ) vL=1 »08 H. ».LL = 1 

P.P I.nT_15 

NPAG = No'aGF(I) • 

.. ...l..= I + 1 . ..... 

GOTO (61 »62»63,64,6F.,66.»67,68).*NPAG 

62. CONTI NUF 

point 102»KPAG 

PR.I N.L..3.2 

print 33 ... 

C .. PRINT PAGF B (PNGLISH) 

_PRINT 21 » ( (MESS (LL) .TIMA (LUf{ ARRAY ( I,LL » 2 ) » I = 1 .08 ) ) » LL=1 »L ) 



I 





PR^I NT 19 

NPATi =_NOAGP(JJ : 

~,i „r L+_-Vl_ ._- :: ^ LLI 

: GO TO (61 »62*6?*6<i,6^»66*67«68 ) «NPA 6 ^ 

- 63 CONTINUE 

‘ ■_ PRINT IO'-^/KPAG ^ 

print 36 ] ^ 1 

PRINt 2 A 

^ print 25 _ 

C ^ PRINT PAGP c (f’ngLTSH) ^ ' _„11 

^ PRINT 23t ( (ME SS(LL) »TI MA(LL) « (ARR A~Y ( I »tf »3 ) >I = l *TO) )_» LL='^L ) 

* print 'i^ ^ ' 

' N PAG = NPA GE fl ) ^ _2 

r=_i_+_i, 

GO TO { 6 1_» 6 2_» 6 3 » 6_5 ♦ 6 6_»_6 7j^6 SJ >NPAG 

,_6A CONTINUE U 

_ print IOAjKRAG 

^ P R I N T" '3 7 

print 34 

print 35 __ . 

i _'c print PAGE^'D { ENGLISH) '* 

~ ■ PR I>^T 2 2 ♦ ( ( ME SS (~LL V* T IM A fCL ) T( A RRAY ( 'ftLl^ATVl = r* 10) ) » LL = 1 ^ Y 

print 19 _ ^ ; '' ' 

' JjNPAG = NPAGE (I ) ^ 1 

i' '=' f +~r 

GO TO (61 »62 ♦63»64,6f',6_6»_67^68 )_»NPA^ -- 

65 continue' ' ^ 1_ r 

^ print' 103 ,'<OAG 1 ■ _ 

,_PR 1 N T 2 fi ^ 

P R I N T' 42 J ^ 

43 _ ' ^ 

C ‘'pr'iNT pTgE F (FnGLISH)_ 

print 26, ( {MESS(LL)_»TIMA(LLI»'(ARRAY( i fLl ♦'5r,I = l,04U »LL=1»L) 
PRINT 19 ~ ' 

NPAG = NPAGE(I)' ' " 

^ - - -- ---- 

GO TO ( 61 ,62 ,63 ,64,6S66,6‘^,68 ) .“NPAG 

* 66"c6ntTnUE‘ 

PRINT l06j<PAG 

, J. ■ ' print' 38 ' ' ■ ■ 


'( 



print 45 

PRINT PAGli F (FMr.LTSH) _ __ . _ 1 • 

PRINT 27, ( ( MESS (J.L1* I IMA aLJ ♦ (A^RAYI I_»LL ,6) ,1 = 1,03) ) ,LL = 1»L) 
print 19 : 

67 CONTINUE ^ 

68 CONTINUE 

"go T0'(1Q9,152_,152)_,NSYS 

152 CONTINUE ] ' ' 

I =_1 _ J 

_ NPAG = NP_AG_EiI) 

!__r= I + 1 ” ; 

G 0_ TO { 8J ,J 2_,^ 7 ,84,8^,8 6,87,88 ) ,NPA6 

8’l’ CONTINUE ' 21 

PRIM T 101 i'K P ^ ^ 

print 30 21_ ^ ] 

PRINT"51‘ 2 " 2 

_ " PRINT PAGF' a ‘(McTr'iC)' _ I 

print _2 0 , '( ( M'ESS( LL ) .TI^MATuL ) , (ARREY( I ,LL , 1 )~, 1 = 1 ,08') ' ) ,LL =Y,L V 
"print 19 ■ 2 1 I 

npag'= NPAGEd'),. " 2.. 2 2 l._l ! 

T ..=. .,1. + 1.. _ 

GO TO .( 8 1., 8 2 , 83, 8 4 ,,.8 6 , 86_,.87_, 8 8 )_j NPA G 

. 8 2 CONTINUE. 

P^.T N T . 10? ,KPAG 2 

PR1NT2 32 ...2 

.PRINT 53. . . ..... 2 

PRINT PAGE 6 (METRIC) 

PRINT 21 , , ( ( MESS ( LDjT IMA (Ll) iTARRFY( I ,LL., 2 ).* I = 1 »08 ) ) , LL = 1 »L ) 

PRINT .1,9 2 ! 

^NPAG_= J^PAGP(n 

,_r2=2 1 +. 2i'”'~’.. _ 

GO TO (’81 ,82,33,84,85,86,87,0'8) ,NPA(i_ 

83 CONTINUE 

PRINT .103 ,<PAG 

PRINT 36.... 2 “ . _ 2.2 

pRiNT_24 2-22 I 22 _ 2___2 

print 56. ... . 

PRINT PAGE C (MPTRIC) 

PRINT 23, ( (MESS(LL) ,TIMA(LL) , (ARREY( I ,LL ,3 ) , 1 = 1 ,10) ) ,LL = 1»’L) 



^ 

'[NPAr, = NPAGF.(.Ii L 

I =. I +.1 J ^ 

. GO TO (81»82»8?»8^,,85»®6.»8X»P-8),»NPA,G 

8^ CONTINUE 

PPINT 10a,<PAG 

P PJ N T 37 ^ 

PPINT, 3 A : - 

print 55 . _ ^ 

PRINT PAGF D (MFfP|C) 

PR INT 22 »_ ( {.MESS { LL ) * T IMA ( L L ) » ( ARRFY ( f*LL » A ) ^ I = 1,10 ).) , LL = 1_,L f 

PRINT. 19_ 

N PAG = N PAGE ( I ) 

I. = .1 „+ 1 

GO TO (81 ,82 ,83, 8 A, 85, 86,67,88) , NPA G 

8 5 .CONTINUE 

PRINT 1.05, <?AG : 

PRINT .2.8 ^ 

.Pp I.N t_A 2 

print 93 _ 

PRINT PAGE E (METRIC) 

PR I NT 26,' ( (MESS'(LL ) ,~T IMA ( LL ) , ( ARREY ( I , U *5 ).,.!.= lVo.AVL»ALf^ 
PRINT 19 --- 

NPAG = _NPAGP(I). ." 1 

L.=_I_+-JL__ : ^ 

GO TO ( 81 ,82.,.83,8A,85 ,P6_,8r,A8_l»M.PAG_ 

R6_C0NTINUE 

print .1.0.6jKPAG 

P.oinT.38 I ^ 

.PRINT. AA 

PRI..NT_.0 5 

C PRINT PAGE F (METRIC) _ 

PRINT 27, ( (MESS(LL)_,TIMA(LL ) , ( A^EY( I ,LL ,6) ,1 = 1,08) ) ,LL=J,L) 

print 19 

87 CONTINUE 'J 

" _ 88 CONT INU'e’_ “ 

19P_C0NTJNUJ 

DO '75“j = l,30 
__ _MFSS(J) = BL<K 
CONT INUE__ 

I PR I NT = ~0' ’ " ■ 

LiNF“ri5__ 

200 'return ' 

END 






i 


l..CRUNGP_ ' SUR, R’uNGF_--- 4TH ORDER .R-K JNTEGRAnON 

SUPPOUTINF RUNRCU rCFI > ^ 

OOU 8 LP PPFCISION. sf : 

DOURLF PPFCISION Do. 

C0MM0N/ST0R8/0MEGA»RF,R.A,PB5 603iH»ALIM L 

COMMON/STCR9/S»MK ^ 

COMMON / STOP 10 / s.t : 

c OMMON / S T OR 1 9 / T 1 M TO ■ ^ ' 

C0MM0N/'st6r25/ LOW , FPT T NY ♦ EPBJ..0 >PT I. ».P I ♦ P. 1 2 ».PT ,H_I GH. 

DIMFNSION S,(9»6)»ST(2.0»1.0.)»K('6.,4) 

DIMENSION DP( 6 )___, ^ 

real. LOW, K 

cl BEGIN RUNGE-KUTTA_INTpGPATIO.M ^ 

.N.^ 1 . j ^ 

IZr_1 ... 3 5 .Q _c' 0 NT I NUE . HI - - - - 

_ timeo = time : 

_ DO 400 J=l ,6 1 

K(J,1) = DT*S( J+'3»1J 

SU,2) = S ( J,1) 4 <(J , l)/2, 

400 CONTINUE..,, . 

_ TIMF * Ti’mE0_-L DT/2. ^ : 

NK = 2 ^ 

.. .. C.ALL ACCEL : 

DO. 4 O 4 J = l ,6 : 

KJ J S.(. j+.3,J?J 

S(J,3) = S(J,1) +.. K( J,2)./2* : 

4 O. 4 . CONTINUE 

time . T I M E P....+_D.T/_2 • 

IZ. 1 ' .. Ic ALL AccE.ril ~~ - Z ^"."""'ZlHZHZHrZZl 1H1._ 

D0..4 0.e_ J..= l.».6 ! 

K IJ.,3 ) , s. DT»S14+3 , 3.) •- 

su,4) = s.( j,i.)_.+_x.( j,3j : 

4 O 8 CONTINUE : 

.J.IMF. .=._LIME0_+_DJ 

_ NK = 4 _ 

CALL ACCE L ' 

60 4ll j = 1,6 _ 

K(J,4) = DT»S(J+3,4) 1_ ” _ ' 

_ _ _ STjJ+10,1 ) = (K(J,1)1+ 2.<K(J»2) + 2.»K(J,3) + K(jt4)_)/6*_ 



4l1L-CDN.T.INU£ 

c STORF, the. SOLUTION ^ 

...DO 490 L=1j6 , 

ST(L»I )_r_.ST(L».Nl_-f_SJJL+TO».I \ 

490 CONTINUE 

C ^ ^ 

<L ^ 

4.9 4_. .1 F,( L . . EO . . 4J .^G.O_T.a_6PO 

I F (.1 . EO , ,_2 I_GO_ TO_49..^ 

GO .TO. 49.7 

.495 CONTI HU E_ 

TIME = TIM E - D T 

_„JpT_=..DI/2a ' 

I..= ^ 

N_=.,l 

GO .TO ,350 

.497.. CONTINUE 

i...=.j,. +_] : • 



00.498 J._.=_.l»6 

S( J.l) = ST(J,3).. .. 

_ 498 CONTJNUE 

NK..=.l 

CALL ACC^.L 

_G.Q_T.O . .3.50 

fOO CONTINUE 1 

DT. .= 2.*DT 

DO 602 J = 4.»6 

.IF(0ABS{<^T(J,2)) .LT. l.OD-7 . ^ 

1 .OR. DAPS(ST(,J»4)1 ,LT. 1 . 0Dr7 ). .GO . T0_.6Q1 

DP(J) = DABSd. O - ST( J .4>/GT( J,2n 

GO TO 6 O 2 

601 .DP( J) .= _Q_,p ^ 

.602 CONTI’^'UE I 

..F.RAC = D^AX] tOP(4)_,nP(6.) ,0P(6) 1 

C _ COMPARE THE DT SOLUTION WITH THE DT/2 SOLUTION 

LELt F.RAC ... LT._E.P_I.mYJ_ilQ_.LQ_6.lJ3 

IF(FRAC .GT, . EPBIG). GO. TO .620. 

_ GO. TO 630 



o.c 


(LID T_ . G T H J r.HI HT = HT r.H 

_GO TO_63q 

62 0 . 1 F ( nT_ .'if. . l6wj_ Gojro _63'0 

n^'£.. =-.Lll^.E_r_i5I 

DT..= DT72. 



I- “ 2 

D0..622.. J = 1 ♦.<? 

S ( J ,1 L . =_,STJ_J ,_1X 

i,2 2. CONTINUE 

DO 621. L=l»6 

S T ( L ♦ 2 ) =. ST.a_._3 1 

6 2 1_ CD NT I N U E. 

GO TO 350 

630. CON.T.I. N'UE ^ 

_ _ IF(DT ,LT. LOW) DT _^LOW 

I .= .2 : 

N_= i _ ; 

STnO.N) = TIME 

. DO 635 M=1 .6 _ 

C ST0RE_THE final SOlUTION_ 

S’T(M.N) = ST(M.4)__+ jLSi(Mj4)_: 

635 CONTINUE _ 

DO 650 ,J = 1.,6 

S (J.N) = s'T(J>N) 

650 CONT.INUE. . 

time. = ST.( l.OjN) 

_ NK = 1 

C_ COMPUTE NFXT AC'cELPFAT ION__ 

CALL.ACCFL ' 

__rptu.rn 

END _ 


- ST(M.2) )/15. 




t 

♦ ♦ 



.CSfTSA*. SUR, SFTSA -^t-VsTOR'f OUTP.UJ/’tPM^'oRA I 

S.UBROUT INF. SEJJ>A^J. 

COr^MON/STORl/ALT,TFf^P,PRF,S»nFNS.VlSC»SO.UND 

. COMwON/STCP2/X,Y»Z»R 

• COMMO'^/STORA/SAREA,MCH»CDD»f:D»MACM»NCri»COS*nRAG : 

l_ __ COMMON/STOR5/TWT»NwT,mdr,v;F1GHT,MASS 

.. COMMON/STOR7/XDOT ,YnOT ♦ZDOT* XDDOT ♦YDDOT.»2DDOJ 1 

COMMON/ STORl 2/ YL ♦ VY.L ip XY ♦GAYL »RLDOT 

.COMM6N/sfOR13/.XL.»2.L ♦.THETA .VXll ♦VZL *ALPH .GAMA.»C.AZ_, SAZ 

COMMON/ STOR 1 5 /.O ♦ VX » VY , VZ ».VT , E.L ♦ P.ANGF.i LA.Tj l6n.G.,.A.CC tGO LAT 

.COMMON/STOR16/SA ♦ARRAY^NN 

COMMON/ STOR 19/ T I M.F, TO 

COMMON/STOR20/ARE/ltIOPTjT.IM, nth, THS. THRUST : 

£_.._.COM.MON/ STOR27yjIPLA.T_, URLQMj.11.PR ♦ 1IRT.I.M 

.COMMON/STOR28/.PLCH , AZ.LCH , ELLCH ♦.VELUCH ,FLCH ,MLCH,ZI.CH,FD LCH, 

1 NDLCH,ZDLCH..„ . 

C OMMON / S T OR 29 / R R A D , A IR AD-J. E L B ADj Vf L.B A D ,.E.R A 0., NRA D. , Z R A O.f E.OP A D i 

1 NDRAD,ZDRAI5_ ^ 

. COM.MON/STOR30/DELTA ,VEL,ACC1 

J? EAL.. 1 1 PLA.T . IJ..PJ-PN , 1Ldp.,.I1PU.M 

.real NRAn»NLCH,NDRAn,Mr)LCH ^ 

. .. RFAL LAMnAO,VA$S , LOW, LONG, LAT, MaCH,MCH 

^ DIMENSION. T ! M { 5 0 ) , T HS ( F 0 ) , T WT ( .SO ) , WPR (.f 0 ) ♦.MCH.( .? 0 Lf.C.DD. (.50 L 

P I me N S 1 ON. S A a 0 , 2., 1.0 1 , A R P A Y..|.10 j.3.0. ♦ 1,01 ^ 

REAL -MASS : 




SA(lfNtJ)..=. ALPH 

SA(2^N,J). = theta 

SA(3tN,J) =. EL..._ 

SA(4,N,J)..= A .Z 

SA( 5tN, J)_ A. ALJI 

SA(6tN,j) = RANGE/6076. 1_155 

SA(.7,N,J). = GDLAT 

SA(R,N,J)_=_LONG 

J. » 2 . _ __ 

SA n ,N,J) = THRUST 

S A ( 2^, N ♦ J ) _E_ W F I GHT 

SA(3fN,J) = DRAG ~ 

SA(4,N,J) = MACH 12 1 _ _12' 
SA(5tN,J) = 0 __ 12 2 * II- 



I 




i 


J) 

Ji 


__SA.t6.»M 
SA(7»M 
SA(fl»N 

- J -Js-3 

SAH sM 
SA(?»N 

_,SA(-3-»N 
.SA(4»N 

- -SA{5»N 
SA(6»N 
SA17>N 

- SA(8»N 

...SAM.O,N»J) 
..-J =.4._, 

SA(1>N 
SA(2»N 
SA(3»N 
..SAi4»N 
S A ( 5 » N 
..SA(6»N 
. SA(7»N 
..SA(8»M 
..SA{0»N 

J = 5 

SAtl »N 
SA{2»N 
SA(3*M 
SA(4»N 

5A(1 »N 
SA(2»N 
SA(3»N 
RETURN 
END 


J.L-J 

JL. 
JL^-MAS*!.. 


.-_ACC_ 

= VT 


= XU. 

=_YL... 


JJ_=_2U 


= -V.YI 

-=-VZl 


J.)..-=..-VXU 

Jl_= 

J). 

JU= RX.Y 

J) = .6AM1 

jj_.=_yEUj:H. 
RLDOI. 


J) = 
J). = 
J) .= 
J.l._= 
J) .= 
J) = 
J) = 
J) s 

J) = 


RRAD/6076, 115.5. 

AZRAO. 

.ELRAD. 

..VELRAD.. 

ER AD/ 60,76..] 1.5.5_ 
NRAD/6076,1155 
2RAD/6075, 115.5. 
EDRAD 
..NDR.AD. 

= ZO PAD 


J) 

J) 

J) 

J) 


IIPLAT 

IIPLON 

IIPR/6076.1155 
IIPTJH 


J) 

J) 

J). 


= VFI 

=..ACC1 


\ 

''•r\ 



I 





:.c. ■ 


4 


% ; 


CTaBLP . SUB. TABLE --- THRUST.WEIGHTtDRAG INPUT. 

SUBROUTINE, table 

__ COMMON/STOR4/SAPFA.MCH»C06»rr>»N'AC^»NCf>.<:r)S*nRAG 


_ COMMON / STOP 5 / TV.'T ♦ N'WT , wpR * WE I GHT 

COMM0N/ST0R6/P0»WRM.WOL»WPO2 »WTM 
COMMON/STOR20/AREA. lOPT. TIM, NTH »THS. THRUST 

common/stor2i/mess,ph.as_. 

COMMON / ST Oi^ 2 4/PHT,MI0 , NTAR,N P H,<P~H 

i 0 FORMAT ( i x\ A6*') I II 

OIMPNSION MeSSOO)_ 



dimension PHT(1 A) ,Nlr)(14) ,WTM(6) ,TSEP(6) »WTP(6) _ 

dimension TIM( 50) ,THS( 50) ,TWT ( 50 ) »WPR ( 50 ) *MCH ( 50J »CDD(5 

REAL LAMOA, MASS, K, LOW, LONG, LAT, MACH ,MCH, LAMDAQ 1_ 

. N AMc L I S T_/ N A M 2./ A P ^ A , T T M , NTH , T HS , I OPT 

’‘NAMPL TSt /NUM'sTfwT ,NWT ,\V0^I_ 

1 N A M P L T S T / N A M A / S_A R F A_, MCH , NC_D , C DO I 

IF(NPH .GT, NTAB") GOlfo" ^0_p_’ ‘ 

NNNN r^NlO.UPH) 

GO TO ' ( 3 A 3’, 3 A A , 3 A 5_,_3 A? J ♦ N NN N I 


?A 3 READ( 5 , 10 ) PHAS 

GO TO 3 AO 

1 -aAA _ CONT I NUE . ' 

C READ IN DRAG TAALF 

RFAO(. 5 , 10 ) _PHAS 

READ ( 5 ,NAMA) HI ' 

GO TO 3 AQ 3 

3 A 5 _C 0 NTINUE _ __ 

C,_ _ READ i’n THRUST r WEIGHT, AND* DRAG Ta'BLE’S 

READ 15 , 10 ) PHAS 

READ( 5 ,NAM 2 ) 

DO 3^6 JslfNTH ' ~ I 

TIM(J) ± TIM( j') P HT(NPH) 

_aa 6 continOe" in 

RPADt 5 ,NAM 3 )' ’ I II 

. . DO 3 A 7 J»l ,NWf" _* _ ,_H_I II HI 

TWT(J) = TWT(j') PHf(NP*H) 

AA 7 CONTINUE ... ' IIH -I_II .. .. I..I 

_WPP? =_ WPP 2 ~ WPP ( 1 ) I IIHI II_ 

I RFAD( 5 InaMA~) 

“ GO TO 3 AR I 

AAO CONTU 4 UE . _IIIII.. 

.rOO return * _ 

END 




( 


Oi 


i 


. CTaBI sup, TABl^ t^PAG COPFFIC lENT . ~ ” . 

. SU.poOUTlNE lABJ 

COMMO^/STOR4/SA9FA,MCH»COD»CDtMACH.NCn.COS»DRAG 

DIMENSION! MCH(50) tf0O(50) 

; . . REAL N'ACH»MCH _ 

- IF (MACH .LE. r,0 T0.116 J 

IF(mach .ge, MCH(\'cr>)j acL-To li? 

DO. ,lpq.js.2.,(aCD 

IF (MACH .LT._ MCHUl )...GO..TO_ll.a 

lOO CONTINUE-^.. 

IlO FRAC a (MACH - MCH( J-1 n/(HCH( J1 - MCH(J-l),) 

. CD a CDD(J-1J + FRAC*(CDD(J)_-_CDD( JrlJ.J 

60 TO 200 

U 5 ..CD_ a^.CDD ( M CD ! 

. 60 TO 200 

116 CD a.CDDd) 

. . ^00 CDS a CD#SAREA 

RETURN 



V 


_C T A B2 . . ” SUr'. . T AB2 J WPU y I N T FRPOL AT I ON^ 

SURPOUT1mF_T.aR? _ __ 

C OMMOM / S T OR 1 / AL T ♦ T F MP ♦ PR 5 S ♦ 0 ENS V V*I s'c ♦ SOUN D _ 

COMMON/STOR6/PC»WRM,woL»WPP2 tWIM ‘ , 

COMVON/STORlR/TIV|r,TO .. .. 

C0MM0N/STCP20/AnEA,10OT,TIM,NTH»THS»THRUST _ 

DIMENSION TTW(50) ♦THS('501 ♦WJV(6) 1 

’ _ IF .( TJ. MC_ ^ LJ . TTMn n r.o TO 1 15 

IF(flMF_,GF, tlM(^THLl_GO TQ 115 

.00.100 j=2»NTH . J 1_ 

..IFfTiMF .,LT..TIM{.J.)i_<10 TO_UQ I 

IpO CONTIN'JF 

.110 TFRAC_s (T'IMF - liM(J-l“))/(TTM(J)-i IMU-l)) ' 

~ ^THPPF = THS(J-i) -I- T FPAC»rTHS ( J ) - THR(J- ln 

GO TO ■(2]0,211 ) . lOpf . 

2lO thrust = THPFF 4 APFA*(PO - PRF'sT _ _1_ 

60 TO poo 

”_.?ii thpust = jhrff «" apfa«~pr'fs' ~~ _H 

GO TO 900 ^ 

1 1 5 thpust j^o , q_ ^ ' 

pOO' RETURN 

FNO 



SUB. TAB3 WPIBHT. interpolation 

..subroutine, tabs : 

COMMON / STOP 5 / TWT ,^M.lT♦••tPR * WEIGHT 

COMMON/STOR6/PO. ;m,wpL>WPP2 »W.TM 

CON'mON/STORIR/TIMP.TO • _ 

O.INPNSIOM TWT(50) «WPR(50) .WTHCe.) 

. I F (.TIME , 6E • TWT < NWT )J. .GO...TOL,.i 1 5 ^ 1 

■DO lOO J=2tMWT 

. I F ( T 1 ME . LT . .TWT IJ ) ). _CO..TO_ 1 1.0. 

CONTINUE... 

IFRAC = (TIME - TWT ( J-1) ) /(TWT(J) - TWT(J-1 
WPP. = WPR( J-l.).'f..JF.R.AC*(WPR{.J)_ - WPRJJ-1))._ 

GO- TO 116 

WPP = WPR (NWT) 

WEIGHT »..WP.P. + .WP.D2_+ .WR_M_..+._WPL. 

..return 

END : 


CTDOUT 


SUP. TPOUT SP^CI&L TAPP OUTPUT 


S'JpeOUTIMP TPOtlT 


COVMOn/STOR16/SA*AoRAv,nN _ 

COMMON /STOR 1 8/ T FP AC » L I NP ♦ P J I ME t CONV ♦ I p ND * l PR T NT *NPAr5E f KP AG 

.COMMON/STOR?A/<LIN 

DIMFNSION TPA(75) »£A(,lO» 2»10) »ARRAY ( 10»30.^10) ^NPAG.E.( 10) L 

...TP A . L I ) =_ PT1ME~ " " ~~ ~~~~~ ~ 


TPA( 
TPA( 
TPA( 
TPA ( 
TPA( 


2 ) 

3) 

4) 
3) 
6) 


= APRAY( 2,L»21._.._. 

= array 


= AopAY(3»L*6J 

= ARRAY (.1jL»6A 


= ARRAY(5,.L*1) 

TPA (..JJ„=_ARRAY1S 


P) 

.R) 


TPA( 
TPA( 
TPA(IO) 
TPAfl?) 
TPA ( 14.) 
TPAdp) 
fPA(20j 
TPA(23) 
TPA (3") 
TPA (55 ) 
TPA{56) 


= AROaY(4,L»2.) 

= ARRAY(. 5,L.,.21. 

= AR0AY(7,L»1) 

= ARRAY( R,L »!.). 

= . array (1 »L»1.) 

„= ..ARR ^‘.Yj. 3 iL * 1 ) 


= array( 4,L rll _ 
= ARPAY( 2»L til,. 
= ARRAY(7,L»2).._ 

= AROAY( 

= ARPAY(2,L»3) 

TP.^i.t.??,) = A RPAY( 3.) .*» ) 

TPA(5P) = appaY(4,L»3) 

ARDAY( 5»L 


TPA(50) = 
TPa(60) = 
TPA(61) s 
TPA( 62) = 
.TpAt.63.J._=.„ 
TPA (64) = 
TPA(65) = 
TPA(66) = 
TPa( 67) = 
fpA(6ft) = 


*3) 

APPAY( 6,L»3) 

ARRAY( 7,L ♦?). . 

APRAY( 8,L,»3 ) 

A R R AXL5.» U.AJ 

ARRAY( 10,L»?1 

ARRAY(4»L»5) 

APPAY(1»L»5) 

AORAY( ?,L»5) 

ARRAY( 3»L»5) 


WRJ TE ( 1 1 )_ ITPAj.JJ . J =1 .71 ) 

DO lOO J=1 ,73 
TPA(J) = 0,0 
lOO CONTIMUF 

RFTURN 

END 




CfoANSr* SU^VItRANSC--- ’in TO POtATING SYS*, 

_____ .SUeP0UTlNF.J?ANSl ( rXfeY^P7./y!Xj_yyj,MZJ I 

DOUBL‘D PPPCI5ION T'^AT^SMAT _ 

C0MM0N/ST0R3/TMAT,SMAT,SPHI /CPHI »CLfSU 

D I mfNS I ON TMAT 13 » 3 ) » SMAT ( 3 ♦ ?J 

WX = TMAT(lfl)*EX + T^•AT(1♦2)*EY + TMATn»3)*f;Z 

WY . = , TMAT(2.»1 )*EX +_ TmAT ( 2 *3) *EY + JMAT ( 2 ».31*E.Z_i 

WZ_ = TNAT( 3«1)*EX + TMAT(3 *2 >* FY -I- TMAT(3»3>»EZ 

RETURN 

^ENT> 



._CTrANS 2 SU3. TRANS2. — ^ ROTATING TO rNEPTIAL Sri* 

SUB.RQUII NE_TRANS.2.( BX *B.yj32LjUXj.yx»itL‘ 


DOUBLE PRECISION TrATjS^'AT 

C0NM0N/STCR3/TMAT,SMAT,SPHI »CPHI fCL»SL 

DIMENSION TMAT(3»3).SmAT(3*3) _ 

UX = SMAT(1»1)*3X + SvAT(1»2)*0Y + SMATn»3)*BZ 

UY = SM.AT(2.»l)f.BX +_SMAT(2»2)*BY + SMAT(2»3)*BZ 

U2 = SM AT<3iT )»BX '+ S^'AT I 3 * 2 ) *Br S MAT(3 >3 I*B Z 

RETURN. : 

END 
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